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Abstract The large amounts of dust detected in sub-millimeter galaxies and quasars at high 
redshift pose a challenge to galaxy formation models and theories of cosmic dust forma- 
tion. At z > 6 only stars of relatively high mass (> 3 Mq) are sufficiently short-lived to be 
potential stellar sources of dust. This review is devoted to identifying and quantifying the 
most important stellar channels of rapid dust formation. We ascertain the dust production ef- 
ficiency of stars in the mass range 3^0 Mq using both observed and theoretical dust yields 
of evolved massive stars and supemovae (SNe) and provide analytical expressions for the 
dust production efficiencies in various scenarios. We also address the strong sensitivity of 
the total dust productivity to the initial mass function. From simple considerations, we find 
that, in the early Universe, high-mass (> 3 Mq) asymptotic giant branch stars can only be 
dominant dust producers if SNe generate < 3 x 10"^ Mq of dust whereas SNe prevail if they 
are more efficient. We address the challenges in inferring dust masses and star-formation 
rates from observations of high-redshift galaxies. We conclude that significant SN dust pro- 
duction at high redshift is likely required to reproduce current dust mass estimates, possibly 
coupled with rapid dust grain growth in the interstellar medium. 
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1 Introduction 



The origin of the significant amounts of dust found in high-z galaxies and quasars (QSOs) 
remains elusive. 

The detection of therr nal dust emi s sion from h i gh-z QSOs at sub-millimeter and mil- 
limeter wavelengths (e.g. JOmont et alll200lll2003l : rCarilli et alll2001bl : iBertoldi and Coxl 



2002 ) indicates far-infrared luminosities > 10 Lq, implying dust m asses of > 10 M, 



and star-formation rates up to 3000 Mq yr ' fe.g.. lBertoldi et a][ 20031). Obs er vational ev 



idence for d ust in these systems has been reporte d bv. e.g.. ^ w ^„ .^^^x., . i ^ 
1 19941) :lLedo ux et al (2002); Priddey et al (2003); R obson et all ( |2004 ); [Ctw v et al (2005 
Beelen et all ( 120061) ; lHines et all ( |2006|) . 

The age of the Universe at z > 6 was less than ~ 1 Gyr. Early star formation is believed to 



Pel et al 



Pettini et all 



have taken place at redshift 10-50 (Tegmark e t al 1997 ; Greif and Brominl , l200^ ; the high- 
est redshift QSO known is at z = 7.1 dMortlock et a]l[201l|) while the earliest observation- 
ally galaxies detected so far include a spectroscopical ly confirmed z = 8.2 gamma-ray burst 



host galaxy dTanvir et~a i [200^:1 ISalvaterra et all l2009l) . a galaxy reported to be at a spectro- 
scopic redshift of z = 8.6 dLehnert et al . 2010h . a gamma-ray burst host ga laxy at z ~ 9.4 



I Cucchiara et all l201lh . and a galaxy at a photometric redshift of z ~ 10 I Bouwens et all 
I2OI l|). The epoch of reionization is determined at z = 10.4 +1.2 dKomatsu et aj 120*101) 



(~500 Myr after the Big Bang). These facts imply that the maximum time available to build 
up large dust masses is at most ~400-500 Myr, and possibly much less. 

Hence, a fast and efficient dust production mechanism is needed. Core collapse super- 
novae (CCSNe) are contemplated to be the most likely sources of dust at this epoch (e.g., 
DwekL 1 19981 ; iTielensl . 1 19981 ; lEdmundsl . I2OOII ; iMorgan and Edmundsl . l2003l ; iMaiolino et al 



20041) due to their short lifetimes and large production of metals. Consequently, several 



theoretical models for dust formation in CCSNe have been developed, which result in dust 
masses of up to 1 Mq per SN within the first ~600 days after the explosion (e.g. . iKozasa et al 



Il989l. 



199ll ; IClavton et aill999l . [200ll ; lTodini and Ferraral.l200ll ; lNozawa et ail2003l) . lDwek et all 



j2007l) argued that 1 Mq of dust per SN is necessary if SNe only are to account for the in- 
ferred amounts of dust in high-z QSOs. 

However, observations of dust in the ejecta of nearby SNe a few hundred days pas t ex- 
plosio n have revealed only ~ 10 ~'*-10~^ Mq of hot (-400-900 K) d ust (e.g., Wooden e Tal 
ll993l ; lElmhamdi et all20()3l ; ISugerman et all2006l ; lKotak et all2009() . Larger amounts (-10"^ 
Mq up to ~ 1 Mq) of cold and warm (20-150 K) dust hav e been rep orted in SNe and SN rem- 
nants (SNRs), afew 10-1000 years after explosion (e.g..[Rho et all2 008. 20091; lOunne etal 
l2009l ; lGomez et all2009l ; lBarlow et all2010l ; lMatsuura et all2oTlh . 

The discrepancy between observationally and theore tically determined dust yields has 
provoked a reconsideration of SN dust formation theories (fcherchneff and D weS 2010l) and 
models incl uding dust destruction have been developed (e.g., Bianchi and Schneideij, I2OO7I ; 
iNozawa et al . 2007; Nath et al, 2008; Silvia et al, 2010). These models demonstrate that dust 
grains can be effectively destroyed in a reverse shock on timescales up to ~ 10"* years after 
the SN explosion. However, they are unable to explain the low observed dust masses at 
earlier epochs. 

D ust production in SNe seems to depend on SN Type (e.g.. lKozasa et a l l2009l ; lNozawa et al 
l2010h . Moreover, intermediate and high-mass asymptotic giant branch (AGB) stars with 



3 



masses between 3-8 Mp have sufficien t ly short lifetimes of a few 10^-10^ years (e.g., 
ISchaller et ail 1991 ISchaerer et all [19931 : ICharbonnel et alL ll993':'Raite ri et alLll996h to be 
potential contributors to dust production in high-z galaxies (e.g., Marche nkd . l2006h . 

In addition to the possible influence from different types of stars on the total amount of 
dust in high-z systems, the prevailing initial mass function (IMF) plays an impo rtant role. In 
the lo cal Universe, an IMF favouring lower mass stars is well established (e.g.. lElmegreenl , 
l2009h while the IMF in the early Universe and in starburst galaxies may be biased towards 
high-mass stars (e.g.. | poane and Mathewsl. fl993l : IPava. ■2008, : ,Dabringhausen et aiL l2009l : 



stars (e.g.,|p 
I et alll201(£ 



iHabergham ( 

In this review we summarize current knowledge about the most important channels for 
stellar sources to produce dust towards the ends of their lives and identify the relevant stellar 
mass ranges contributing to the total amount of dust in galaxies. We determine the ranges 
of dust production efficiencies of AGB stars and SNe and address the influence of various 
IMFs on the dust productivity of stars between 3^0 Mq. Based on this insight we review 
what is currently known about the stellar contribution to dust in high-z galaxies. The review 
is arranged as follows: we first summarize our knowledge about the late stages of stellar 
evolution of massive stars (Sect.[2]l. In Sect. [3] some fundamentals of dust grain formation 
and characteristics are described. We address the complexity of determining the amount of 
dust theoretically and observationally in Sect. |4] (evolved massive stars) and Sect. [3] (SNe). 
Dust production efficiencies are quantified in Sect.|6]and the impact of the IMF on the total 
dust productivity is discussed in Sect.[7l The inference of large amounts of dust in massive 
high-z galaxies and QSOs, along with theoretical models addressing this topic, is reviewed 
in Sect. [8] Sect.|9]provides a summary of the main conclusions of this review and an outlook 
for future directions. 



2 The late stages of massive stellar evolution 

For the most likely dust producers, such as AGB stars and CCSNe, the majority of the 
dust production takes place at the end stages of their evolution. Therefore, pertaining to the 
observed presence of dust in galaxies and QSOs at z > 6, only stars which live short enough 
to die before the age of the Universe at this redshift are conceivable sources of dust. 

In Fig.[T]we illustrate the relation between the minimum zero-age main sequence (ZAMS) 
mass of stars and the redshift at which they die. We have considered three different epochs 
for the onset of star formation. For a formation redshift of z = 10 we find that the lowest 
mass of a star to be a potential source of dust at z = 6 is 3 Mq. Less massive stars can be 
excluded because their lifetimes are longer than the age of the Universe at this redshift. The 
effect of the metallicity with which a star is bom is small. 

Owing to this ascertainment, we are solely interested in the high-mass (> 3 Mq) stellar 
population. We therefore briefly summarize what is known about the end stages of massive 
stellar evolution, which eventually govern the dust production of these stars. 



2.1 The first stars 

The first generation of stars, so-called Population III (Pop III) stars, played an important role 
in reionizing the Universe and were responsible for the early enrichment with metals. They 
are belie ved to have formed i n dark-matter mini halos of ~I0^~* Mq at redshift z ~ 10- 
50 (e.g.. iTegmark et all [l997l : lO'Shea and NormarJ. |2007|) . The very first stars (Pop III.l) 
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Fig. 1 Relation between stellar mass, stellar lifetime and redshift in the early Universe. The graphs show 
the minimum ZAMS mass of a star dying at a given redshift for an onset of star formation at three different 
epochs: z = 15 (dotted curves), z = 12 (dashed curves) and z = 10 (solid curves). The colour coding 
corresponds to different metallicities: Z = 0.001 (black) and Z = 0.040 (magenta). The vertical dashed line 
marks a star dying at z = 6, similar to the highest-redshift QSOs known. The grey sh aded region indicate s 
stars with masses be twee n 2-3 Mp. The metallici ty dependent lifetimes are taken from lSchaller et ai il992h . 
ISchaerer et a] <1993h and lCharbonnel et al il99A The cosmological model used is a ACDM Universe with 
f/o = 70 km s"' Mpc-\Q^ = 0.3 and = 0.7. 



formed in isolation and are expected to have been relatively rare, only about 10% by m ass of 
all generations of Pop III stars (e.g. JGreif and BrommL[200^ : lMcKee and TanLl2008h. From 
simulations it is predicted that these stars are very mass ive, ~ 10^"^ Mp (e.g., lAbel et 
|2002; Bromm and Larson, 2004; Schneider et al, 2006; Y oshida et all 20061) . The formation 
of the second generation of stars (Pop III. 2) is influenced by the radiative and mechanical 
feedback ef fects of t he first stars and is found to be delayed by about ~ 200 Myr (e.g., 
I Johnson et al. 2007; Yoshida et ^l2007ah. The critica l mass of t hese s tars is suggested to be 
lower, about ~ 30^0 Mm (e.g., 'Yo shida et alLl2007bl ; lNormanl . l2010l) . For a review on the 
first stars we refer the reader to iBromm et a]| l2009h . 

According to lHeger et all (l2003h . stars with metallicity Z = and masses between 40- 
140 Mq and above 260 Mq collapse into black holes, while stars in the mass range 140-260 
Mq die as pair instability SNe (PISN e). The explosion wi l l entir ely disrupt the star, leaving 
a quite peculiar chemical signature jHeger and WooslevL [20o3) which is manifested in a 
strong odd-even effect of the produc ed nuclei. So far, on ly one supernova, the Type Ic SN 
2007bi has been reported as a PISN lOal -Yam et all 200ft) . However, theoretical modeling 
indicates that SN 2007bi may also be consistent w ith an energetic core-collapse SNe with a 
main sequence progenitor mass of ~ 100-280 Mq (iMoriva et ail2O10l ; I Yoshida and Umedal . 
|m1 1). The typical PISN signature ex pected to be observable in the first and most metal-poor 
stars has not been d etected yet (e.g., [ Beers and Christliebl[2003) . The reason for the non- 
detection is unclear. lEkstrom et all ( 20081) discuss the possibility that, under the conditions 
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of fast rotation, CNO line driven Wolf-R ayet (WR) win ds and magnetic fields, the PISN 
stage of very massive stars can be avoided. iKarlsson et all ([2008) argue that stars formed out 
of gas enriched by primordial PlSNe are more metal-rich and thus the signature of PlSNe 
would not be expected in metal-poor stars. 

The metal enrichment by Pop 111 stars leads t o formation of low-mass Pop 11 stars, as 
soon as a critical metalli city of ~ 10~*-10~^ dBromm and Loebl . l2003l : [Schneider et all 
l2006l : lTumli^^ . l2006h is reached. This transition is expected to take place fast due to a 
rapid meta l enrichmen t ( Maio et al |2010|) . although metal-free regions survive over longer 
timescales. lGreifetall(l2010h showed in a cosmological simulation that one single PISN (at 
z - 30) can enrich the mini halo in which it forms uniformly up to Z = 10~^ Zq and induce 
Pop 11 star formation (at z - 10). 

The IMF for the first stars is con sidered to have been very top heavy with high char- 
acteristic masses > 35-100 Mq (e.g.. lBromm et a 1 12OO2I : iTumhnsonl . l2006l : lYoshida et all 
|m08). For the generation of Pop II stars, top heavy IM Fs with som e what lower character- 
istic masses or Salpeter-like IMFs are usually assumed. iTumlinsonI ( |2006|) points out that 
metal-free star formation is relatively scarce at redshift z ~ 6. Owing to the above discus- 
sion, PlSNe and the very first stars are disputable to be major dust sources for dust-rich 
galaxies at z ~ 6. 



2.2 AGB stars 

Stars in the AGB phase are in their late stages of evolution. They have initial masses in the 
range ~ 0.85-8 Mq and have completed the helium-burning phase in their centers. AGB stars 
have low surface temperatures (max 3500 K) but high luminosities (a few times 10^ Lq) and 
have built up so-called helium- and hydrogen-burning shells around their degenerate cores 
of carbon and oxygen. The hydrogen burning shells deliver the energy needed to maintain 
the high luminosities. During the AGB evolution the stars develop quite strong winds with 
incre asing mass-loss rates toward s their late stages whereby they lose some of their matter 
(e.g.. lSch6ier and Olofssorll200lh . The very late stages are characterized by intense mass- 
loss, which increases towards the end, when th e stars enter a super-wind phase with mass- 
loss rates up to 10""* Mq yr"' (e.g. jBowen and Willson. 1991; Schoier and Olofsson, 200l[). 
In general, the stars lose up to ~ 80 % of their masses during the AGB phase and form 
circumstellar envelopes of gas and dust. Low and intermediate mass stars (at the lower mass 
end of the AGB mass range) end their lives as white dwarfs. However, the final fate of stars 
with masses around 8 Mq might be different (see Sect. 12.4b . 

AGB stars can be broadly divided into three distinct classes based on low-resolution 
spectra: (i) The oxygen-rich M-stars whose spectra are dominated by bands due to TiO 
molecules, (ii) the carbon-rich C-stars whose spectra are dominated by bands due to C2 
and CN molecules and (iii) S-stars which are ne ither rich in oxygen nor ca rbon, identified 
through their strong bands due to primarily ZrO dLattanzio and Woodl . l2003l) . 

The distribution of C- and M-stars is a function of stellar mass and initial metal abun- 
dance. For low initial metallicities it is eas ier to form C-stars as l e ss car bon needs to be 
dredged up. Stellar evolutionary models by iKarakas and Lattanzid ( |2007|) predict that for 
Large Magellanic Cloud (LMC)-like metallicities, M-stars evolve from low (1.0-1.5 Mq) 
and high (5.0-8.0 Mq) mass stars, while C-star s originate from intermediate (1. 5-5.0 Mq) 
mass stars. The latter has also been found by e.g.. lVassiliadis and Wood d 19931) and lZiilstra et all 
1 I2OO6I) . 
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Generally, AGB stars with masses above 4 Mq experience hot-bottom burning (e.g., 
iBlocker and Schonbemed . Il99ll : IP'Antona and MazzitelliL Il996h leading to a reduction of 
the a mount of carbon which can b e dredged up. Models of massive low-metallicity AGB 
stars I Ventura and D'Antonal , l2009h show that Z = 10"'* stars appear as C-stars during most 
of their AGB phase while a slightly higher metallicity of Z = 6 x lO"'* yields M-stars. 



2.3 Core collapse supemovae 



CCSNe are divided into two different classes. Type I and Type II, and their subtypes (e.g., 
iFilippenko. 1997, and references therein). Type II SNe are defined by the presence of hy- 
drogen lines in the optical spectra while Type I SNe are defined through their absence. The 
CCSNe subtypes can be aligned roughl y in the order of in creasing progenitor mass, starting 
with II-P, II-L, Iln, lib, lb and Ic (e.g.. [Anderson and Ja mes. 2008). As we discuss below, 
there is no one-to-one correspondence between progenitor mass and spectral Type. However, 
the alignment of the SN subtypes might correspond to increasing mass-loss of the hydrogen 
envelope (for Type II SNe) and subsequent stripping of the he lium envelope (Type I SNe) 
of the progenitors (e.g.. lNomoto et allll995l ; lMaund et aill201 1 *). The main characteristics of 
the subtypes leading to the typical CCSN classificat ion scheme are sum marized in Table[T] 
The most common CCSNe are Type II-P SNe. ISmartt et all ( |2009|) find, that the mass 
range of the progenitors of Type II-Ps is betwee n 8.5"^]'° and 17 ±1.5 Mq. A lower mass limit 
of ~8 Mq was found independently by And erson and Jamea 1 2008]'). However, theoretical 
predictions from stella r evolution models (e.g., Heger et alll2003l : lEldridge and Toutil2004l : 
IPoelarends e"taill2008h indicate a higher upper mass limit for II-P SNe: For solar metallicity 
it is ~25 Mq and increases with decreasing metallicity. It is therefore unclear what happens 
with stars more ma ssive than 17 Mp. The progenitors of II-Ps are found to be red super- 



giants (RSGs) (e.g.. IWooslev and Weavel 1 19861: IS martt', '2009*: 'C rockett et all 1201 ih . Con- 

rtt et aT 



firmed e xamples include SN 2003gd jSmartt et all l2004: Maund a nd Smarttll2009l) and SN 
2008bk jMattila et alLl2008bl . 1201(1 : IVan Dvk et aiBoiOli. The Type II - P SN 2002hh likel y 
arose from a RSG progenitor of around 16-18 M^ ( IPozzo et all l2006l : ISmartt et ai. l2009h . 
However, R SGs up to 25 M© have been detected in the Local Group (e .g.. iLevesgue et al 
2005, 2006). One pos sibility i s that they collapse and form a black hole JSmartt et alll2009l : 



Heger et al, 2003; F rver et ai 12007 ). In that case, they either appear as very faint SNe or 
no explosion is observed at all due to fallback of ^*Ni. Alternatively, more massive stars 
might end their lives as Type II-L SNe. A possible exa mple is SN 2009kr whose progenitor 
could be a yellow super giant (YSG) of mass 15-24 Mq jElias- Rosa et al, 2010; Fraser e tal, 

Imid). 

Stars more mas sive than 25 Mq evolve into WR stars and likely explode as lb or Ic SNe 
jMassev and Olseil 2003: Cro wtheiL l2007h . A small fraction of these involve a relativistic 
explosion ( Sode rberg et alLl2O10h leading t o broadlined Ic SNe, sometimes accompanied by 
a gamma-ray burst (e.g.. lHiorth et a 112003). During their precursor luminous blue variable 
(LBV) stage, they lose their hydrogen envelope either through massive eruptions or in peri- 
ods of enhanced mass-loss and form a rather de nse circumstellar disc (C SD). The WR phase 
lasts for approximately 10^ years (e.g.. lMevnet'an d Maeder, 2003; E ldridge and Vinkl . l2006h 
whereupon the star finally explodes as a CCSN, leaving a black hole. 

Stars more massive than 25-30 Mq may explode during their LBV phase before entering 
the WR stage. T his was the case for S N 2005gj (Kotak and Vink, 2006; Trundle et al, 2008) 
and SN 2005gl iGal-Yam et alll2007h. which b oth appeared as very bright CCSNe of Type 
Iln. The Type Iln SN 2005ip jSmith et all2009l) had a different progenitor, probably a RSG 
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Table 1 Definitions of different types of core collapse SNe^ 



SN Type 


Defining 
characteristics 


Progenitor 
mass range 
in Mq 


Progenitor 
characteristics 


Prominent 
examples 


Type 11 


hydrogen present 








il-r (plateau ) 

TT T i]inf^r\^ 

ii-L t^nncdij 

Iln (narrow line) 

Iln 
Dn 

nb 


blue, almost featureless 

spectrum 

very blue, almost 

featureless spectrum 

narrow emission lines on a 

broad base 

low luminosity 

very luminous 

similar to Type lb SN, 

little hydrogen present 


(7) 8—25 

~ 1 S 9S 

— LD — ZJ 

~ 8-10 

> 25-30 

> 25-30 


RSG 
1 oyj 

SAGB 
LBV 

WR, binary 


oJN IVoyL, 3W ZUUigCl 

CM OflfiQlrr 

SN 2008S 

SN 2005gj, SN 2005gl 
SN 1993J, SN 2008ax 


Type 1 


hydrogen deficient 








lb 
Ic 


helium rich 

helium deficient, no Si 11 


> 25 

> 25 


WR 

WR, binary 


SN 1999dn 
SN 2007gr 



Note s. ^The typ es of CCSNe are classified based on their spectral appearance and ph oto metric evolutio n 
(e.g.. iFilippenkd 1 1 997V Information on the progenitor is taken from lSmartt et al i2009ll and lSmarti i2009l) . 
^The lightcurve exhibits a 'plateau' for an extended periode after maximum brightness. "^The lightcurve (in 
magnitudes) declines linearly with time. 



of roughly 20^0 Mq. This shows that Type Iln SNe may arise from either stars with LBV- 
like mass ejection if they are very luminous, or from massive RSGs with a strong wind 
interaction if of moderate luminosity. 

The fate of stars more massive than roughly 17-25 Mq is sensitive to mass-loss ef- 
fects, depending on magnetic fields, metallicity, binarity, or rotation, although the de tails of 
these dependencies are not well understood (e.g., see for a review |PulsetalLl200i) . As a 
consequence, there is no simple relation between SN Type and progenitor mass. It seems, 
however, that these stars rather explode as I ln, lib, lb or Ic SNe tha n ordinary Type II-P SN. 
The progenitor of the Type lib SN 2008ax dCrockett et aj . 2008; Pa storello et ail2008h was 
a late-type 28 Mq WR star with strong nitrogen emission lines in the spectra (a so-called 
WN L star). The mass of the progenitor of the Type lib SN 2003bg was estimated to be 20-25 
Mo llMazzalieta]L[200^. The Type lb SN 1999dn seems to be consistent with a progenitor 
of mass 23-25 Mr. /Senet ti et allbouh . An example of a Type Ic SN is SN 2004gt, where 
the progenitor mass is estimated to be > 40 Mq (Maund et al, 2005). For the Ic SN 2002ap 
a single star pr ogenitor of 30-40 Mq has been proposed, but with very high mass-loss rates 
llCrockettetalLl2007h . 

Binary interaction between two lower mass stars has been contemplated as the progen- 
itor for some Type lib and Ic SNe. For the lib SN 1993 J a companion star was clearly 
observed ((Maund et al, 2 0041). Other exam ples where th is scenario has been invoked in- 
clude the Hb SN 20 lig ( iRvder et alLl2006l). SN 2q02ap jCrockett et alll2007t) , SN 2004gt 



dMaund et a]ll2005l) and SN 2008ax 



Crockett etaill2008h . 
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2.4 The critical progenitor mass range 8-10 Mq 

The fate of stars in the mass range ~ 8-10 Mq is ambiguous since the mass border between 
high-mass AGB stars and CCSNe is smeared out. Moreover, the decisive facto rs for the 
development of stars in this mass range are unfortunately rather uncertain (e.g. jNomotd . 
[T984, 1987). We denote this range the critical progenitor mass range. 

Stars in the critical progenitor mass range may either evolve directly into II-P SNe or 
form an electron-degenerate core of oxygen, neon and magnesium (0-Ne-Mg) and enter the 
super AGB (SAGE) phase. SAGE stars can either become 0-Ne-Mg white dwarfs or turn 
into electron capture SNe (ECSNe). 

The appearance of an ECSN could be as a faint and ^^Ni poor II-P SN such as SN 1994N, 
SN 1999eu or SN 2005cs jPastorello et"aiLl2004l . [200^ . However. ISmartt et all (l2009h did not 
find any signature or convincing evidence that faint and ^^Ni poor II-P SNe are ECSNe. The 
inferred luminosities for progenitors in this mass range rather favour normal Type II -P SNe. 
A peculiar example is the faint and ^^Ni poor II-P SN 2007 od i Andrews et ail2O10h . whose 



spectra favour a SAGE star as the progenitor llnserra et alll201 1 ). 



Alternat i vely, a ECSN could occur as a low-luminosity Type Iln SN such as SN 2008S 
jPrieto et ai l2008h . The progenitor mass of SN 2008S was deterrnined to be -6-10 Mp 
and the progenitor could have been a SAGE star jPrieto et a i l2008l : lEotticella et all [20091 : 
IWesson etajbOlOh . 

According to lWanaio et"al ( l2009h . about 30 % of all CCSNe co uld appear as ECSNe if 
all stars in the mass range of 8-10 Mq end in the explosion channel. IPoelarends etal ( l2008l. 
and references therein) suggested that onl y the most massive (9-9.25 Mq) stars will explode 
as ECSNe, representing ~ 4% of all SNe. isiessI ( I2007l . l2008h showed that at very low metal- 
licity (Z = 10"^) , the m ass range for stars becoming an ECSN is much broader (7.6-9.8 M©). 
[Thompson etaj il2009h suggested the rate of ECSNe to be ~ 20% of all CCSNe. 



2.5 Type la supernovae 

Type la SNe are characterized by the absence of hydrogen and presence of silicon in their 
spectra and are be l ieved to result from a thermonuclear exp losion of a carbon-oxygen white 
dwarf fe.g.. lLivid.l2000l : lHillebrandt and Niemeveil.l2000h . An explosion takes place when 
the white dwarf reaches the Chandrasekhar mass through external mass supply. However 
the nature of the progenitors and explosion patterns are controversial. Two scenarios are 
currently favoured; (i) a single-degenerate model where a main sequence or giant compan- 
ion star transfers mass by Roche lobe overflow dWhelan and IbeJ. Il973l : Ipihk et all l2007h 
and (ii) a double degenerate m odel where the companion star is also a white dw arf and the 
two objects merge (I ben and Tutukovl , ll984l : [Webbinkl . ll984l ; iPakmo r et al, 2010). The mass 
range of stars possibly exploding as Type la SN is 3-8 Mq (e.g., iMaoz . .2008) which means 
that the stars become C-0 white dwarfs after having evolved as AGE stars. This gives rise 
to rather long delay times between the formation of the progenitor sys tem and the explosion 
due to the long lifetimes of these stars. From explo sion models (e.g.. lGreggio and Renzinil 
ll983l : lMatteucci and Recchill200ll : rGreggioll2005h different delay times are predicted. Ob- 
servations indicate that there may be two different progenitor channels, resulting in SNe 
la with delays of either < 40 Myr or > 2.4 Gyr since progenitor formation (Erandt etal 
l2010l) . lMannucci et all l l2006h suggest that half the SNe la explode already after about 100 
Myr while the other half have longer delay times of about 3 Gyr. 
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Table 2 Main characteristics of some prominent dust species 



Dust species 


Chemical 


[K]" 


Spectral characteristics. 


Reference 




definition 




prominent bands [pm]^ 




Carbon-rich environment 


Amorphous carbons, a-C:H 


sp^/sp', H 


> 1700 


^ 0.2, 3.4, 6.85, 7.25 


1-8 


PAH' 


fusions of C(,li(,^ 


< 1700 


0.2-0.26, 2-50 


7, 9, 10 


Graphitic carbon 




~ 1600 


-0.22 


11, 12 


Silicon carbide 


SiC 


> 1700 


- 10-13 


12-14 


Oxygen-rich environment 


Olivine 


[Mg, Fe]2Si04'" 


~ 1300 


-0.7-1.5,9, 10-11.6, 18, 20 


15-19 


Forsterite 


Mg2Si04 


~ 1300 


- 10, 11.3, 69 


15, 17, 20, 21 


Fayalite 


Fe2Si04 


~ 1000 


- 10.6, 11.4, 93-94, 110 


15, 17, 20, 22, 23 


Pyroxene 


[MgFe]Si03^ 


~ 1300 


- 10-20, 40.5 


15,24 


Enstatite 


MgSiOj 


~ 1300 


- 9.7, 19.5, 26-30 


15,24 


Magnetite 


Fe304 


~ 800 


- 17, 25 


15, 25 


Corandum 


AI2O3 


~ 1700 


broad at - 13(12.5-14) 


15, 26 


Spinel 


MgAl204 


~ 1200 


- 0.3, 0.5, 2, 13, 17, 32 


15, 19, 27 


Calcite 


CaCOs 


~ 800 


-6.8, 11.4,44, 92 


28 


Dolomite 


CaMg[C03]2 


~ 800 


-6.6, 11.3,60-62 


28 


Iron 


Fe 


-900 


featureless 


15, 29, 30 



References. (1) 'Jager et aj fmsd); (2) 'Duley and Williams' (1981'); (3) 'Jones and NutW ('2011'); (4) 
Dartois et al (2004); (5) Frenklach et al ( 1989); (6) Pascoli and Polleux (2000); (7) Jager et al (2009a); (8) 
Mennella et al ( 1997); (9) Cherchnefi (201 1); (10) Cherchneff et al 1 1991); (1 1) Draine and Lee ( 1984); (12) 
Sharp and Wasserbura ( 1995); (13) Lodders and Feglev ( 1995); (14) Mutschke et al ( 1999); (15) Gail ( 2010); 
(16) Koike etal (1993); (17) Jager etal (1998a); (18) Pitman et al (2010); (19) Zeidler et al (2011); (20) 
[Koike et al (2003); (21) Koike et al (2010); (22) Hofmeister ( 1997); (23) Pitman and Hofmeistei (200 6); (24) 
iKoike et all feOOOh ; (25|Koike et a] (1981); (26 Koike et all 119951) ; (27) iFabian et all I i2001i) ; (28) iPosch et all 
<2007l) ; (29) |Palikl <1985l) ; (30) Kemper et al (2002); 

Notes. "Tc is an approximate dust condensation temperature which varies with different properties, such 
as the proximity to thermodynamic equilibrium, the prevailing pressure or the composition of the gas (e.g., 
ISharp and Wasserbursill99^ ; lGaill . l2010l) . ''The exact peaks depend in most cases on properties of the grains 
such as temperature, crystallinity, size, morphology or impurities. Most silicates and oxides ai'e highly trans- 
parent in the ultraviolet (UV), visual and NIR and absorption data are either hardly available or laboratory 
data are not reliable since the values are often too small to be measured with standard methods of spectroscopy 
jZeidler et alll201 lb . "^Polycyclic aromatic hydrocarbon (PAH), "^benzene (C^Hf;) is an aromatic hydrocarbon 
(AH) and PAHs consist of several fused aromatic benzene rings, e.g., pyrene (CieHio), pentacene (C22H14), 
coronene (C24H12) ^Olivine and pyroxene ai'e non-stoichiometric compounds with the chemical formulation 
for olivine: Mg2j:Fe2(i-j^)Si04 and for pyroxene: MgvFe(i_v)Si03, with < x < 1. 



Having discussed the end-stages of massive stellar evolution, we next address the dust 
formation processes associated with massive stars (Sect.[3l|4]and[5]l. 



3 Fundamentals of dust formation and grain species 

Dust is formed by a series of chemical reactions in which atoms or molecules from the 
gas phase combine into clusters of increasing size. The molecular composition of the gas 
phase determines which atoms and molecules are available for grain formati on and grain 
grow t h. The sizes of dust gra i ns are in the range of a few 0.01-1 jim (e.g., iMathis etal 
ll977l : IWeingartner and Drainell200lh . 
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The dust formation process is typically described as a two-step process, i.e., the con- 
densation of critical seed clusters out of the gas phase and the subsequent growth to macro- 
scopic dust grains of certain sizes and species. The nucleation proce ss in the ma jority of 
current models is based on the so-called classical nucleation theory which 
was developed to explain the formation of water droplets in the Earth's atmosphere. It has 
been found that at temp eratures between ~ 70 K and ~ 2000 K and d ensities in the range 
~ 10"'^-10"'^ e cm~^ jFedeii 1 1964 IClavtonl [197^ : ISedlmayj, 1 19941) thermodynamically 
stable clusters can form. Note, however, that the applicability o f this theory to astrophysical 
environments has been questioned (e.g.. lDonn and NiithLll985h . An alternative theory based 
on chemical kinetics has also been applied to describe dust formation in diverse environ- 
ments (e.g., Frenklach and Feigelson, 1989; Cherchneff et al, 1992). 

Dust grain formation depends on various critical parameters such as for example the 
sticking probability, a. This parameter depends on, e.g., the material under consideration, 
the internal energy of the grains, the impa ct energy or the temperature of the gas. However , 
the exact st icking probability is uncertain ( lDraineLri97^ : lLeitch -Devlin and Williamsl[T98^ ; 
lGai ll2003l) and is therefore often taken to be unity, for simplicity. In this case all colliding 
particles stick together, leading to a maximum amount of dust to be formed under the given 
nucleation and growth conditions. 

The resulting dust species depend on the environment. Stellar environments are mainly 
rich in either carbon or oxygen and depending on the most abundant element, predominantly 
either carbonaceous dust or silicates will form (see Table|2]for an overview of some common 
dust species). 

Carbonaceous dust mainly consists of the element carbon (C) and has manifold forms of 
appearance defined by the types of carbon hybrid orbitals (sp, sp^ and sp^) leading to differ- 
ent bond structures between the C atoms (e.g., Henning et al, 2004). Amorphous carbon is 
typically characterized by the ratio of sp^ to sp^ hybridized bonds. Graphitic grains are com- 
posed solely of sp^ and nano-diamonds of sp^ hybridized bonds, respectively. Dependent on 
the ratio of sp^ to sp^ and the impurity of other elements in the single bonds (i.e., H, N), var- 
ious subtypes of either amorphous c arbon, such as hydrogenated ttmorphous carbon (a-CiH) 



or nano-diamonds are created (e.g.,lDulev anc 



iHenning etaill2004 l Jones and D'Hendecourd 



Will 



ariisl 1 1 98 ll : iMolster and Watersl 12003 



20041) . In an environment which is not only 



carbon-rich but also H-rich, polycyclic aromatic hydrocarbon (PAH) can form at low tem- 
peratures (Tc < 1700 K) (e.g., Cherchneff et al, 1991; Jiiger et al, 2009a; C herchneff. IQV^ 
Typic ally, PAHs are fusions of several aromatic benzene (CgHg) rings (e.g.. lCherchneff et all 
Il992h and constitute the building blocks for the condensing solid particles, so-called soot 
grains (e.g., Jager et al, 201 1). The largest PAH molecul es condens ing together with the soot 
grains are found to have 222 C atoms (C222H42) (e.g.. jjager et all l2009ai bL and references 
therein). 

Other dust grains found in a carbon-rich environment are for exampl e FeS, MgS or 
different p olytypes of sili con carbide (SiC), such as ff-SiC or /?-SiC (e.g., iBorghesi et all 
T985.;,OrofinqetalL ll99lh. which condense at high temperatures (e.g., > 1700 K) (e.g., 
Daulton et all 2002 , and references therein). 



Silicates are the most stable condensates. Typically, silicate grains consist of [Si04]'* 
tetra hedra in conjuncti on with Mg^"^ or Fe-"^ cations. For reviews on cosmic silicates we re- 
fer to lHennind l l2010alBb . Overall, silicate grains can be categorized into amorphous silicates 
or crystalline siUcates. The crystalline lattice structure allows the tetrahedra to share their 
oxygen atoms with other tetr ahedra. This leads to the formation of different types of silicates 
jMolster and Kempe^l2005h , such as forsterite, fayalite, olivine, enstatite, ferrosilite or py- 
roxene (see Table [2] for details). Forsterite and enstatite are the most abundant crystalline 
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silicates jMolster and WatersLl2003l . and references theirein). Carbonates are chemical com- 
pounds of the characteristic carbonate ion C0^~ with elements such as e.g., Ca, Mg or Fe. 
Common examples are calcite (CaCOs), magnesite (MgCOs), dolomite (CaMg[C03]2) or 
siderite FeCOs . Other dust grains formed in an oxygen-rich environment include corundum 
(AI2O3), grains of the group of spinels (i.e., spinel (MgAl204), magnetite (Fe304)), silica 
(SiOa) or metallic iron (Fe). 

4 Dust from evolved massive stars 

4. 1 Dust from AGB stars 

In the local Uni verse, AGB s t ars are the prime sources of dust injected into th e interstellar 
medium (ISM) llGehrzLll989l : ISedlmavilll994l : lDorschner and Hennind.ll995l) . The dust is 
injected as part of the intense mass-loss during the late stages of AGB stellar evolution 
(see Sect. 12. 2t . The driving mechanism of the mass-loss is believed to be a combination of 
ther mal pulsation and radiation pressure on du s t grain s resulting in slow dust driven winds 
re.g.. lH6fner e "tIiLll998l:lHii fner and Ander i^ l2007h with typical velocities between 3-30 
km s~'. 

The dust composition in AGB stars depends on the C/0 ratio in the photosphere of the 
star which is directly connected to the nucleosynthesis in the stellar interior. Newly formed 
elements like carbon and oxygen are mixed to the surface by a deep convective zone. The 
mixing processes o ccur during the thermall y pulsing AGB (TPAGB) phase and also involves 
the external layers dlben and Renzinilll98ll) . The TPAGB ph ase lasts approx imately lO'*"^ yr 
depending on the stellar mass and number of thermal pulses ( lB16ckeilll995h . The stellar pul- 
sations cause atmospheric shock waves propagating through the atmosphere. Subsequently, 
gas is lifted above the stellar surface, producing dense, cool layers favourable for possible 
solid particle formation (e.g.. lH6fner et allll998h . The ongoing nuclear burning and dredge- 
up changes the relative abundance of carbon and oxygen as the stars evolve. A change in the 
C/0 ratio results in a change of the spectral type (see Sect. I2.2| | and the composition of the 
dust (see Sect.[3j. 

- M-type: C/0 < 1 results in an oxygen excess since all carbon is bound in CO molecules 
creating an oxygen rich environment where either silicates or carbonates are formed 

- S-type: C/0 » 1 leads to an exhaustion of C and O which are almost completely bound 
in CO. For this type no abundant grain forming elements are available and grain species 
are defined by the less abundant elements. 

- C-type: C/0 > 1 creates a carbon rich environment (all oxygen bound in CO) where 
predominantly hydrocarbon molecules and carbonaceous dust forms together with some 
silicon carbide. 

Deriving the dust driven mass-loss characteristics of AGB stars is difficult and the 
current understanding is based on numerical models. Detailed time-dependent dynamical 
models featuring a frequency-dependent treatment of th e radiative tr ansfer have success- 
fully explained the mass-loss mechanism for C-stars (e.g.. lH6fner et a i|2003;Hofner,2doi; 
[Winters et a ilaooO). In such stars, amorphous carbon grains form from the excess of carbon 
at high temperatures. Mass-loss is enhanced by the radiation pressure on such grains which 
efficiently accelerates the dust particles away from the star, dragging the gas along. Models 
invol ving C-rich dust driven mass-loss are well tested and ar e consistent with observations 
(e.g.. lGautschv-Loidl et aill2004lNowotnv et alll200ll201(]h . 
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In the cas e of M- s tars the oxygen environment leads to formation of preferentially Fe- 
free silicates dWoitkel . I2OO6I : iHofner and Andersenl . l2007h . in particular olivine and pyrox- 
ene type grains. S uch grains are consis tent with observed features in infrared (IR) spec- 
tra of cool giants jMolster et all I2002al lbll3). However, small (< 1 yum) grains of Fe-free 
silicates would result in insufficient radiation pressure to drive a wind due to thei r trans- 
par ency at wavelengths correspo nding to the flux maximum of AGB stars. iHofnej ( l2008h 
and 'Mattss on and Hofned ( I2OIII ) have shown that larger grains of sizes in a very narrow 
range around 1 /im can drive a wind. This grain size range is also consistent with grain sizes 
observed in the ISM. 

Dust formation and mass-loss in S-type stars pose substantial problems l H6fneij . [2009h . 
According to observations, S-stars show circumstellar physical properties similar to C- and 
M-stars (e.g.. iRamstedt et al l2009h . However, the equality between the most abundant el- 
ements O and C inhibits the formation of known mass-loss driving dust species such as 
amorphous carbon or micron-sized silicates in sufficient abundances. Several minor dust 
species have been proposed which however are either not abundan t enough or of too low 
opacity to enhance mass-loss (e.g.. lFerrarotti and Gaiil2002l . I2OOQ) . While some of these 
species possibly play an important role, it remains uncertain which dust types or processes 
drive the winds. 

AGB stars are important suppliers of molecules w hich play a crucial role in forming 
dust and are important for mass-loss mechanisms (e.g., loiofssorj. Imi . 1 1 997I : Ik^^ et all 
Il998l . and references therein). IR and sub-millimeter observations of AGB stars hav e re- 
vealed many different t ypes of molecules (e.g., CO, SiO, R\Hs, etc.) (e.g., |jus ttanont et all 
ll996l : lYange"tIill2004l). In theoretical stud i es of M-, S- and C-stars , using a chemical ki- 
netic approach (e.g., Oierchneff et~^ Il99ll . 1 19921 : ICherchnefi[2003) . a wide variety of the 
observed molecules could be reproduced. Apart from the molecules formed only in specific 
environments (C- or 0-rich), species such as CO, SiO, HCN and CS have been found theo- 
retically and observationally in all types of AGB stars (e.g., ICherchnefd . [200^ : lDecin et all 
l2008l . l20I(]h . 

Studying the influence of metallicity on the mass-loss and dust formation processes in 
AGB stars is important t o und erstand their role in the early Universe. Theoretical investi- 
gations bv lWachteretdl ( l2008h showed that the wind velocity decreases with lowering the 
metallicity of low-mass AGB stars, while the mass-loss rate remains unaffected. The latter 
also resulted from models by Matt sson et all (200^ under the condition that the amount 
of condensable carbon in low-metallicity AGB stars is comparable to that of th e more 
metal-rich counterparts. Using the James Clerk Maxwell Telescope, iLagadec et 
performed CO observations of six carbon stars in the Galactic Halo and the Sagittarius 
stream and came to similar conclusions: The mass-loss rates of C-stars are unaffected by 
metallicity but the expansion velocities for metal-poor C-stars are lower. Spitzer Space Tele- 
scope observations of the LMC, the Small Magellanic Cloud (SMC) and the Fornax Dwarf 
Spheroidal indicate that mass-loss rates for M-ty pe stars are more sensitive to metallic- 
ity while metal-poo r C-stars are unaffected (e.g ., Ziilstra efa ^ 20061 : iGroenewegen et al 
I2OO7I : iLagadec et al l2007l : iMatsuura et al 120071 : ISloan et all l2009l) . The amount of dust 
produced by M stars is found to decrea se with decreasing metalli city while for C stars 
it remains unchanged ( IGroenewegen et ak ,2007, : ,Sloan et al, 2008(). Ivan Loon et a 1 ( l2008h 
present ESOA^LT spectra of a sample of dusty C stars, M stars and red supergiants in the 
SMC. A comparison of the properties of molecular bands to similar data in the LMC indi- 
cates that dust formation in M-stars as well as in C-stars is less efficient at lower metallicities. 

Typical mass-loss r ates obtained observationally and theoretically are between I0~^ 
and 10"^ Mo yr"' (e.g.. ISchoier and Olofssonl . I2OO1I : IWiUsoJ. I2OO7I : iMattsson et all l2008l : 
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iMatsuura et all l2009h . The mass-loss rates and the dust-mass-loss rates are li nked via the 
eas-to-dust mass ratio, for which a canonical value of 200 is often assumed (e.g. jLagadec et a! 
|2009; Sloan et al, 2009). Estimated values of the gas-to-dust mass ratio obtained from CO 
observations coup led with either radiative tra ns fer or dynamical mod eling range between ~ 
200 and 700 (e.g. jGroenewegen et alLflQQSbl fl lRamstedt et ail2008h . 

It is important to stress that the models discussed above are developed for low and in- 
termediate mass stars. A theoretic al dust formation model fo r AGB stars in the mass range 
of 1-7 Me has been developed bv lFerrarotti and Gaill 1200^ . Dust yields are calculated for 
several metallicities and result in total dust masses up to a few times 10"~ Mq. The con- 
sidered grain species are silicates, iron dust, SiC and carbon, which are the most abundant 
grain types in M-, S- and C-type AGB stars. The model combines synthetic stellar evolution 
models with a non-equilibrium dust formation prescription. The dynamical treatment of the 
stellar outflows is simplified in that stationary flows are assumed and hence the mass-loss 
rate is an in put parameter because it c annot be determined self-consistently. Nevertheless, 
the model of lFerrarotti and Gaiilbooeh is currently the only available source which provides 
dust yields for AGB stars covering a large range of stellar masses and metallicities. 



4.2 Dust from red super giants and Wolf-Rayet stars 

RSGs are evolved O or B stars at the He-burning stage and arise from massive stars (< 
40 Mo). Stellar evolution models bvl Mevnet and Maedeil ( l2003h sugges t that the RSG phase 
lasts for about 0.4 Myr for a ~ 25Mo star or 2 Myr for a ~ IOMq star. lMassev et all ( |200^ 
estimated a dust production rate of about 3 x 10"* M© yr"' kpc"^ for RSGs in the solar 
neighbourhood which is about 1% of the dust return rate of AGB stars. 

For RGSs with a luminosity < 1000 Lq no evidence for dust production has b een found 



in studies of the globular cluster 47 Tuc |Bover et a iL l2010l : lMcDonald et alLl2oTlh (contrary 
to the findings of lOriglia et 1ij200l l2010^ However, stars more luminous than ~ 1-2 xlO 



Lo do appear dusty (e.g..lMcDonald era]ll2009l . [2CiTl[) . 

From observations of the H II region NGC 604 in M33, it was found that RSGs ap- 
pear more extinguished than WR stars, indicating large amounts of dust around the RSGs 
( lEldridge and Relafld.l2oTTI) . WR stars are the successors to massive RSGs which undergo 
strong winds which drive away the dust created prior to the WR phase ( lEldridge et a]Ll2006f) . 
Thus WR stars appear less extinguished. 

Although the above examples show that dust is produced around evolved massive stars, 
it is unclear how much of the dust survives the subsequent SN explosion. There is evidence 
for ongoing dust destruction in the expanding SN bla st wave of the SNR Cas A, as well as 
dust evaporation due to the UV flash from the SN jDwek et ail2008h . 



5 Supernova dust 

As outlined in the following, on average, a few times 10"* Mq of relatively hot dust (~ 500- 
1000 K) has been reported from CCSNe at early epochs while large amounts of cold dust 
(< 50 K) have been claimed in SN 1987A and SNRs which are a few 100-1000 yr old. 
In contrast, theoretical models predict that a high amount of dust in the SN ejec ta can 



form within the first 600-1000 days (Kozasaetal, 1989, 1991; Clayton etal, 19991, 20011 



9^ 



iTodini and Ferraral , l200ll : lNozawa et a]| . i2003> ; iBianchi and Schneideri . i2007i ; iCherclmeff and Dwekl 
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l2010h . The calculated dust masses are of order 10"'-1 Mq for SNe in the mass range 12^0 
Mq, for metallicities between 0-1 Zq. 

Pertaining to this controversy, we next address the difficulties of deriving dust masses 
in SNe from either theory or observations and provide a status of the current observational 
situation. 



5.1 Theory 

The fir st m odels to investigate dust formation i n SNe were developed by e.g.. lCemuschi et all 
d 1965b and lHovle and Wickramasinghd (Il970h . More recent work has addressed the various 
dust species eind amounts of dust formed in SN ejecta, but the models are still in their infancy. 
The applied theories (standard nucleation theory or chemical kinetics) for dust formation 
in SNe is similar to that used in models for ot her stellar environments such as (i) stella r 
outflows of AGB stars (Sect.|4lJ or LBVs (e.g..lFerrarotti and Gai]Ll200lLlGail et ail2005h . 



but also (ii) su bstellar atmospheres (e.g.. iHelling et all l2008h oTTiii) brown dwarfs (e.g., 
lBurrowsl[200^ and references therein). 



In addition to the uncertainties in the applied dust formation theories (Sect.[3]l, the SN 
dust models are hampered by complex physical processes which are not well understood, 
such as the SN explosion and subsequent expansion of the ejecta. The amount of dust and 
the variety of dust species formed in the theoretical models thus strongly depend on the 
assumptions made. 



5.1.1 Dust formation models based on standard nucleation theory 

iTodini and Ferraral ( l200lL hereafter ItFOiI) investigated the formation of dust in Type II SN 
arising from progenitors with 12-35 Mq and metallicities between zero and solar. The nucle- 
ation of dust grains is based on the classical nucleation theory. For the formation of CO and 
SiO molecules chemical equilibrium is assumed. The ejecta are considered to be spherically 
symmetric and the chemical elements are fully mixed. The gas temperature and density are 
uniform throughout the considered volume. The temporal evolution of the temperature is 
defined by the assumption of an adiabatic expansion of the ejecta. For the kinetic energy of 
the explosion two different values are considered. In most models, amorphous carbon grains 
are typically the first grains which condense out of the gas phase about 300-400 days past 
explosion. Large seed clusters made of monomers are able to condense and amorphous 
carbon grains grow to large grain sizes of about 300 A. 

Most of the amorphous carbon dust is formed at a gas temperature in the ejecta of about 
T = 1800 K. As the ejecta expand other dust species condense at lower gas temperatures, 
i.e., corundum at T ~ 1600 K, and then magnetite, enstatite and forsterite at T ~ 1100 K. 
Typically fairly small dust grains of about 10-20 A of these species form. At zero metallicity, 
and in the lower energy case, the calculated total dust masses per SN are in the range 0.08 
Mq < < 0.3 Mq, but are increased when a higher explosion energy is assumed. The dust 
masses per SN increase with increasing metallicity, e.g., the amount is three times higher 
for Zq relative to zero metallicity. The obtained log-normal grain size distribution is found 
to be rather insens itive to metallicity. 

The model of ItFOII has been revisited by iBianchi and Schneidej ll2007l) to study the 
effect of a reverse sh ock on the dust grains. Another grain species, Si02, is added to those 
already considered bv lTFOll Moreover, only clusters with a minimum number of monomers 
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of either N >2 or N > 10, and discrete accretion of these, are considered. These modifica- 
tions lead to an alteration of the log-normal grain size distribution of all dust species except 
for amorphous carbon grains and result in a larger mean grain size (and less numerous 
grains). With increasing A', less Si-bearing grains of large grain sizes form while amorphous 
carbon grains are not affected. In the case of solar and sub-solar metallicity, around 0.1- 
0.6 Mq of dust is formed per SN. However for Z = no dust is produced for progenitors 
more massive than 35 Mq. The final dust masses per SN are sensitive to varying the sticking 
probability a between 1 and 0.1. Assuming a = 0.1 leads to significantly reduced total dust 
masses (0.001-0.1 Mq of dust for progenitor masses below 20 M©). Higher a leads to larger 
grains. Si-bearing grains are significantly more affected by lower values of a than amor- 
phous carbon grains, and in some cases the amount of dust in Si-bearing grains becomes 
negligible. A shift of the size distribu tion of carbonaceo us dust grains towards larger grains 
with higher a has also been found bv lFallest et ailll201lh . 

iBianchi and Schneided ilOOl h used a simple semi-analytical model t o treat the dynamics 
of the reverse shock. The model is based on analytical approximations bv lTmelove and McKed 
lll999h for the velocity and radius of the forward and reverse shocks in the non-radiative 
ejecta-dominated phase and sub sequen t Sedov-Taylor phase of SNRs. For the energy and 
ejecta mass, values similar to the lTFOll formation model are adopted and three different val- 
ues for the ISM density are investigated. Furthermore, a uniform density distribution inside 
the spherically symmetric ejecta is assumed along with a uniform distribution of the dust 
grains. The grain size distribution is considered to be the same throughout the ejecta. It has 
been found that due to erosion caused by thermal and non-thermal sputtering, the grain size 
distribution is shifted towards smaller grains. Depending on the density of the ISM about 
2-20 % of the initially formed dust mass survives (higher fraction at lower density). About 
4-8 xlO** years after explosion the reverse shock has penetrated 95 percent of the original 
volume of the ejecta. 

iNozawa et a hereafter lN03h studied dust formation in the SN ejecta of zero 

metallicity stars (13^0 Mq) taking also PISNe (with masses of 170 or 20 Mq) into ac- 
count. The classical nucleation theory is assumed, but following iGail et all ( Il984l) . a non- 
steady state nucleation rate is calculated. The temporal evol ution of the density and temper- 
ature are calculated following hydrodynamical models bv lShigevama and Nomotol ( Il990l) 
and a multifre quency radiative tran sfer code coupled with the energy deposition from radia- 
tive elements dlwamoto et alL[2000l) . respectively. The models distinguish between unmixed 
ejecta and mixed ejecta. In the unmixed case, the ejecta are divided into different layers, 
each of different elemental composition, i.e., the innermost regime (consisting of elements 
such as Fe, Si and S) followed by oxygen-rich layers (composed of elements such as O, Si 
and Mg), and outermost, a He layer. The mass of each layer varies with progenitor mass. 
In the mixed case, all elements are assumed to be uniformly distributed. For either case, a 
formation efficiency of unity is assumed for the key molecules CO and SiO, and the total 
amount of freshly formed dust is found to increase with progenitor mass. The total amount of 
dust per SN for the mixed ejecta generally is found to be larger than for the unmixed ejecta. 
For SNe between 13^0 Mq about 2-5 % of the progenitor mass condenses into dust while 
the corresponding fraction for for PISNe between 140-260 Mq is 15-30 %. In the mixed 
case the ejecta are oxygen rich due to the assumption that the formation of CO molecules is 
complete. Consequently, only oxide grains such as forsterite, corundum, enstatite, Si02 or 
magnetite condense. The most abundant grain species for SNe are SiOz and forsterite. In the 
unmixed case various different grain species condense in each layer depending on the ele- 
mental composition of those. The main grain types formed are carbon, Fe, Si and forsterite. 
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The average grain radius of each grain species depends on the elemental composition and 
the gas density at the formation site. 

Similar to the study of iBianchi and Schneided ( l2007h . iNozawa et~a il l l2007h investigate 
dust destruction caused by the impact of the reverse shock in the SNR phase of zero metal- 
licity stars. The ejecta are assumed to be spherical and to expand into a uniform ISM with 
primordial composition, where three different cases for the hydrogen number density are 
co nsidered. For the density an d velocity structure of the ejecta the hydrodynamical m odels 
of lUmeda and Nomotd (l2002h together with the dust models (mixed and unmixed) of lN03l 
are adopted. Three diiferent radiative cooling processes are included. Dust destruction by 
sputtering and the deceleration of dust grains due to gas drag are taken into account while 
the effect of charge of the dust grain is neglected. 

Initially, very large grains (> 0.2 yum) are found to be expelled into the ISM through 
the forward shock while their size is only marginally reduced through sputtering. Smaller 
grains are either destroyed through sputtering in the post-shock flow or are trapped and 
remain behind the forward shock. The critical grain size below which dust particles are fully 
destroyed is sensitive to the density of the ISM and is found to be in the range 0.01-0.2 jim 
for a hydrogen number density in the range 0.01-10 cm~^. The grain size distribution of 
the surviving dust is therefore dominated by large grains. The fraction of dust destroyed is 
found to be higher for the mixed grain model than for the unmixed, as the mixed model lacks 
grains larger than > 0.01-0.05 fim. Furthermore, the final fate of the dust grain s depends 
on the thickness of the hydrogen envelope of the progenitor star jNozawa et all [2010). In 
the case of a thin hydrogen envelope (as expected for Type lib SNe) smaller grains form. 
Moreover, the reverse shock encounters the ejecta much earlier than for SNe with a thick 
hydrogen envelope (as is the case for Type II-P SNe). In the latter case the reverse shock 
encounters the dust ~ 10^""* yr after explosion, depending on the density of the ISM. 



5.1.2 Dust formation model based on chemical kinetics 



Models for dust formation in the SN outflow o f a zero metallicity star of 20 Mq an d PISNe 
(170 and 270 Mg) have been accomplished bv lCherchneff and Dwek' ('2 009l.l2010h . For the 
temperature and density structure of the SN ejecta, the models of NOS! are adopted. The 
temporal evolution of those quantities is calculated by assuming that the ejecta follow an 
adiabatic expansion similar to the mode ls described above. T he ejecta velocity is for sim- 
plicity kept constant and a mixed ( Umed ~and Nomotol . l2002l) and unmixed case ( N03) are 
considered. It has been argued that the commonly adopted assumptions of thermodynamic 
equilibrium as well as the standard nucieation theory are ina ppropriate for describing dust 
forma tion in the dynamical flows of SN ejecta (see Sect. [3ll. ICherchneff and Dwekl 1 20091 . 
l2010h therefore use a chemical kinetic approach for the formation of molecules and dust 
grains. The chemical kinetic description of the ejecta is based on (i) the initial chemical 
composition of the gas and (ii) a set of chemical reactions describing the chemical processes 
in the ejecta. This new approach leads to smaller dust masses by a factor of ~ 5 a nd to a 
diffe re nt chemical compositi on of the formed dust compared to the models of either ItTOiI . 
lN03l or lSchneider et a il i2004t) . The most abundant grain species which form in these models 
are pure silicon, silica and silicates, while carbon dust is negligible. 
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5.2 Inferring dust masses from observations 



Deriving the mass of dust from observations is equally complex. Warm dust emits in the 
near-infrared (NIR) and mid-infrared (MIR) wavelength range, whereas the emission from 
cold dust is shifted to far-infrared (FIR) or sub-millimeter wavelengths and is often difficult 
to differentiate from cold foreground material. In addition, it is impossible to infer the struc- 
ture of dust grains and their spatial distribution within the ejecta from observations. Hence, 
the derived dust masses rely on the models and techniques used to interpret the data. 

The methods mainly used to infer the existence of dust are based on observations of 
either (i) the attenuation of the red wings of spectral lines at optical/NIR wavelengths or 
(ii) the thermal emission from dust grains. Most observations of SNe and SNRs have been 
made with the Spitzer Space Telescope since its launch in 2003, because ground-based MIR 
observations are difficult. The instruments onboard of Spitzer cover in the w avelength range 
3.6-160 yum. Earlier observations, e.g., of SN 1987A l lWooden etailll993h were performed 
with the Kuiper Airborne Observatory operating in the 1-500 pm spectral range. The start of 
operation of the Herschel Space Observatory in 2009 has facilitated FIR and sub-millimeter 
observations (detectors sensitive to wavelengths between 55-625 ^m ) of SNe and SNRs 
as has already been accomplished, e.g., for Cas A and SN 1987A ( iBarlow et all l20IOl : 
lMatsuuraetall201lh . 



5.2.7 Dust masses in SN ejecta 



The attenuation of broad and intermediate spectral emission lines, e.g., the He I, Ca II IR 
triplet or O I line, is a relatively reliable and usually pronounced signature of the presence 
of dust. Using this method direct confirmation of ne wly form ed dust in the ej ecta has bee n 
presented for so me SNe, including SN 1987A (e.g.. iDanziger e t al, 1989; Lu cy et a idii), 
SN 19901 (e.g.,|Elmhamdi et al', '20 04)), SN 1999 em (e.g., Immhamdi et ak i2003h . and SN 
2004et (e.g.,|Sahu et al, 2006; Kotak et al l |200^ . Evidence for formation of new dust not 
only in the ejecta but also in the post-shocked shell of Il b/IIn SNe was rev ealed for example 
forSN 1998S (e .g., Pozzo et al 2004.), SN 2005ip (e.g., ISmith etlill2009i). SN 2006ic (e.g 
Smith et all l2008b; Mattila etail2008al) and SN 2007od (I Andrews et alLl2O10l : llnserra et all 



20111) . Unfortunately, it is difficult to quantitatively derive the amount of dust, or its compo- 



sition or geometry (e.g., jCotak, 2008), with this method. 

Thermal emission from dust is typically detected in IR observations of SNe as a NIR or 
MIR 'excess'. Such an 'excess' may arise from newly formed dust in the SN ejecta or in the 
cool, dense shell of post-shocked gas within the forward and reverse shock. The new dust 
may be collisionally heated by hot gas in the reverse shocks, or heated due to radioactivity 
or optical emission from circumstellar interaction. Alternatively, thermal emission could 
be caused by pre-existing dust in the circumstellar medium. In this case the dust is either 
collisionally heated by hot, shocked gas, the flash from the SN or it is heated due to the 
interaction between the ejecta and the circumstellar matter The latter two cases result in an 
'IR echo' due to light travel time effects. It is challenging to differentiate between newly and 
pre-existing dust from observations of thermal emission, although thermal emission caused 
by an echo seems to appear at earlier epochs than emission due to dust formation, which 
takes place a few hundred days past explosion. 

Studies of SN light curves are also useful since, in case of an echo, the light curve shows 
characteristic features. However, either scenario might contribute to the l ate-time IR flux as 
was the case for SN 2004et, SN 2006ic, SN 2007 it and SN 2007od (e.g.. lKotak etaill2009l ; 
iMattila et alll2008al : [Andrews et alll201oll201 Ibl) . 
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5.2.2 Dust masses in SN remnants 

In old SNRs (see Sect. l5.3.T] | it is possible that most of the dust is cold and has escaped de- 
tection in MIR studie s. Sub - millim eter observat ions with SCUBA h ave been accomplished 
for the Cas A ( Dunne et all l2003h and Kepler ( iMorgan et al l2003h SNRs. The first mea- 
surements resulted in very large derived dust masses (~ 0.3-3 Mq) at cool temperatures of 
about 17-18 K. However, in particular for Cas A it has been suggest ed that most of the 
sub-milimeter emission h kely arises from foreground molecular clouds ( IKrause et a 
IWilson and Batrl£j.l2005h . Similar considerations and new calculations led to a downwards 
revision of the dust mass for Kepler of about a factor of two jGomez et all l2009l ). Using 
either sub-m illimeter p ola rime try (Dunne et al, 2009) or the Herschel Space Observatory 
instruments iBarlow et all 12010) . lower dust masses were obtained for Cas A as well, al- 
though for either remnant the obtained amount of dust is well above the average results of 
MIR-studies in SNe at early and late epochs. 



5.2.3 Caveats and uncertainties 

Deriving the dust masses in SN ejecta and remnants is basically similar to the method used 
for deriving the amount of dust in galaxies (see Sect. 18. lb . Based on the method discussed 
by Hildebrand ( 1983), the dust mass for a single temperature component is determined from 
the flux density observed at some frequency, X, as 

F(/l)Z)2 

Mi = , (1) 

Ki(A,a)B(X,Tiy 

where F{X) is the total flux, Z)l is the luminosity distance to the object, is the dust 
temperature, B{X, T^) is the black-body Planck function and Ki{X, a) = (3/4)g(/t, a)/(pa) is 
the dust absorption coefficient for a (spherical) grain type. Here Q(A, a) is the dust absorption 
efficiency, p is the dust bulk density and a is the dust particle radius. The temperature can 
be derived from a spectral fit. The luminosity of a single spherical grain of radius a and 
temperature is given as Ld(^) = A-na^nB{A, Ti)Q{A,a). 

In the Rayleigh limit, a < A, the absorption coefficient k is independent of the particle 
radius a, thus k = Ki(A), which is usually adopted since exact grain sizes and grain size 
distributions are unknown. 

The main uncertainties in deriving the dust mass are (i) the considered dust species, (ii) 
the optical constants, i.e., the dust absorption coefficients for the considered dust species 
and (iii) the unknown grain size distribution. For reported dust mass estimates the adopted 
dust grain composition often varies. In addition, different optical constants are often applied 
for s imilar dust species i.e., ( i) for graphite and silicate grains the o pti cal consttints are taken 
from lDraine and Led ( ll984) ,|Draing ('1985^. Ossenkopf et all (ll992h orLaor an d Drainel ( Il993h . 



(ii) for Mg protosilicates from IDorschner et al tl980a) o r iJager e t al ( 2003) ar id (iii ) for 
amorphous carb on grains values fro m Hannei ( Il988h or iRouleau and Mar tir] 1 1991 ) are 



assumed. While ISouchet et 

3 (2004*) preferred silicate dust for SN 1 987A, Ercolano et all 
^200^,) adopt large amounts of graphite grains. For m ore recent SN e, Fox et al (2010) rule 



out silicates and use only graphite for SN 2005ip as do lMattila et a\ ( l2008ah for SN 2006jc. 
[Andrews et al tOld) and [Andrews et a favour an amorphous carbon dominated 

model for SN 2007it and SN 2007od. Spectroscopic evid ence for sihcate d ust was re- 
vealed through a large, but declining SiO mass in SN 2004et jKotak et"aill2009|). For Cas A. 



iHines et aill l2004l) adopted a magnesium protosilicate-based grain model from lDorschner et all 
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lll980lj) and lArendt etail ( Il999h . iRho et all ( l2008h fi t th e spect ra with a variety of different 
grain species b ased on the t heoretical models o fiNOalandlTFOl' but favour magnesium proto- 
silicates, while lDunne et al ( 2003, 2009) assume grains which are either amorphous or have 
a clumpy, aggregate structure. Silicates and graphite dust have been assumed in the Crab 
and the SNR B0540 tGreen et alLl2004l : lTemim et all 1 20061 : 1 Williams et ail2008h . Assuming 
silicate rather than graphite dust typically leads to higher inferred dust masses. It is evident 
from the above highlighted examples that the determination of the grain type composition 
is not trivial. 

Further complications in deriving the amount of dust from SNe arise from ambiguous 
considerations about the SN ejecta physics. In most cases it is unclear whether the ejecta are 
mixed or unmixed, and additionally a uniform dust and gas distribution is often assumed, 
while there seemingly is evidence for mixing and clumpy ejecta, as we explain below. Mix- 
ing in the ejecta likely can be explained by the theoretically observed instability of the 
nickel bubble during explosion of the SN leading to Rayleigh-Taylor instabilities forming 
in the post-shocked ejecta (e.g., Chevalier and Klein. 1978l: lAmett , ll988l : rHerant and Benzl 
1 19911 : iHerant and Wooslevl . 1 1994 : .Kifonidis et aL ,2003i') . This might also support sugges- 
tions for the presence of undetected larger a mounts of dust at early epochs, if dust pains 
are assembled in opt i cally thick clumps (e .g., Lucv et al, 1989, 1991; Elmhamdi et al, 2003; 



IWooden et al 1 19931: ISueerman et alll2006 : Ercolano et al, 2007 ; Meikle et al, 200?,.) . Ac 



cording to Meikle et all 1 2007h . dust in the ejecta of SNe can become optically thick in the 



MIR for dust masses exceeding a few times 10""' Mq. However, in most of the cases where 
clumpy models have been applied, significantly larger dust masses than for smooth models 
were not found. In particular also these models f ail to explain the large dust masses pre- 
dicted by theo retical models dWooden et"aill993l : lErcolano et aj l2007l : Meikle et ailioOTi : 
I Andrews et ai l2bl0. 201li). 

A scenario of dust grain growth in SN and SNRs over longer timescales of a few 10- 
1000 yr could explain the difference in dust mass at early and late epochs. Once a stable 
cluster has formed, further growth to macroscopic dust grains can take place. The growth 
regime of dust grains extends to lower temperatures and densities than for the nucleation 
regime (see Sect.|3]l. However, significa nt growth is re s tricted by t he ava ilable condensable 
material and dilution of the SN ejecta jPraind . 1 19791 : ISedlmavl Il994h . The amount and 
timescale of grain growth might also be dependent on the Type of the SN, examples of 
which are the SNRs B0540-69.3, Cas A or the Crab nebula. For the latter, an extended dust 
grain growth phase could possibly explain the presence of large dust grains jTemim et all 
l2006h . 



5.3 Quantitative evidence of dust from supernovae 

Bearing in mind the caveats discussed in Sect. 15.21 we now proceed to summarize current 
reported observational evidence for dust arising from SNe. 

5.3.1 SNe and SN remnants with reported dust properties 

Direct evidence for dust formed in SN ejecta and remnants has been reported for only a few 
cases so far. 

- In the peculiar Type II supernova SN 1987A at most a few times 10~^-10~^ Mq of 
dust at epochs between 615-6067 days past explosion was found jOwek et"ai 1 19921 : 
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I Wooden et aiLll993l : lBouchet et all I2OO4I : lErcolano et alLl2007l) . FIR and sub-millimeter 
observations with the Herschel Space Observatory 8476 and 8564 days after its ex- 

i ilosion reveal about 0.4— 0.7 Mp of cold (17-23K) dust formed in the SN outflow 
Matsuura et alll201l1 ; lLakicevic et allbOllh . 
- At epochs between 214—1393 days past explosion, dust masses of at most a few times 
lO"** Mq at temperatures of a few hundred K was inferred for the Type II-P super- 



2007od (e.g., Elmhamdi et al, 2003; Suserman et alL 


2006: Meikle et aL 2007; KotaS 


2008 


; iKotak et all 20091: Andrews et all l20ld. 2011b 


: Meikle et all 


201 ll: Szalaietall 


2011 


). 



X 10 ^ Mo and 2.0 x 10 ^ Mq on day 678 with a clumpy model. In contrast, [Meikle et all 
( l2007h inferred only 4x10' Mq of hot dust and concluded that the mid-IR emission 
from this SN cannot support a dust mass of 2.0 x 10"^ Mq. They also argue that the 
difference in the result s may be due to the pre sence of a larger component of cold dust 
in the smooth model of ISugerman et all ( |2006|) . 
- A quite peculiar case is SN 2006jc. Two years before explosion, a LBV-like outburst 
was detected and ass ociated with the progenitor of SN 2006jc (Nakano et al, '2006"; 



Pastorello et alll2007 ) , whic h has been suggested to be a very massive star CFoley et al, 
2007l ; |Pastorello et alLl2007h . dence for ongoing dust formation in a C SD behind the 
forward shock already at 55 days after explosion was reported by e.g., loi Carlo et all 
(I2OO8); Smith et al (2008b), but just a modest amount of 3 x lO"** Mq of dust was in- 
ferred | Mattila et al, 2008^. Interestingly , also larger dus t masses of ~ 8 x 10~^ Mq 
(iMattila etaLl2008an or ~ 3 x 10"^ Mq (S akon et ail2009h condensed in the mass-loss 
wind of the progenitor prior to explosion was observed. 

In SNRs, at an age of a few 100-1000 yr, larger masses of rather cold dust seem to be 
present. For example, observations of the SNR Cas A result in a few times 10~' Mq of hot 
dust (>170 K) and a few times 10"^ Mq of warm and cold dust (< 150 K) for the entire 



SNR (e.g . lArendt et all [l99^ : jPouvion et all l2001bl : iHines et all l2004l : iKrause et all l2004l : 



aNK (e.g . .lArenc 

iRho et all 12008V An amount of ~ IMq of dust at a temperature of ~ 20 K was recently 
suggested bv Dunne et al (2009i) . Observations with the Herschel Space Observatory result 
in a resolved cool dust component (~ 3 5 K) in the unshoc ked interior of Cas A with an 
estimated mass of 7.5 x 10"^ Mq of dust ('Barl ow et aibOlOh . For the SNR 1E0102.2-7219 
observations by Sandstrom et al (2008) have shown that 3 x 10"^ Mq at 70 K are present 
as newly formed dust in the ejecta, which has already encountered the reverse shock. From 
observa tions in SNRs arising f rom Type II-P SNe such as B0540, SN 1987A or the C rab 
nebula dwilliams et a l. 2008; B ouchet et aiLl2004l : iOreen et al, 2004 : Temim et all l2006h an 
average of a few times 10"^-10~^ Mq of dust has been inferred. 



5.3.2 Type Iln supernovae and LBVs 

There is growing evidence for dust from Iln SNe and LBVs. SNe of Type Iln arise from 
stars at the lower mass end of CCSNe (8-10 Mq) or from stars with higher masses in con- 
nection with LBVs (> 20 Mq). In either case they have undergone strong mass-loss and are 
surrounded by a dense and hydrogen-rich circumstellar disc. 

In the case of ECSNe (see Sect. l2.4l i appearing as Type Iln SNe, dust formation seems to 
be quite efficient. SN 2008S was embedded in a dust enshrouded circumstellar shell and the 
progenitor was likely a SAGB star. The dust enshrouded phase lasted for ~ 10'* years prior 
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to explosion (iThompson etalLl2009h and could be associated with the super-wind phase of 
SAGB stars. 

The Type Iln SN 2005ip is an e xample where dust was formed in the post-shocked shell 
l lFoxeta]ll2009l :[ Smith et a i|200i).A pre-existing large dust shell containing ~ 1-5 x 10 
Mo of warm (~ 400 K) dust and a hot (~ 800 K) dust compone nt of about 5 x 10"* Mq 
arising from newly formed dus t in the ejecta has been found bv IFox etall ll2010h . For the 
Type Iln SN 1998S. |Pozzo et all ( |2004) inferred a dust mass of > 2 x 10"^ Mq. 

SN 2006gy was classified as the most luminous Iln event known dOfek et all [20071 : 
[Smith et al, 2007; Kawabata et al, 2009), but the explosion mechanism remains uncertain 
( iOfeke t al, 2007; S mith et alil2007l,l2008ij,l2010ij) . NIR observations two years past explo- 
sion i Miller et ali . l20fol) showed a growing NIR excess which can be explained by a massive 
shell of around 10 containing around 0.1 Mq of dust heated by the SN. The existence 
of a dusty shell has been proposed to be due to LBV eruptions lasting over ~1500 years 
prior to the SN explosion ( iSmithetall[2008al^ . The large mass of the circumstellar medium 
(CSM) of ~ 1 0-20 Mq and the lik ely SN ejecta mass of 10-20 M© require a progenitor mass 
of ~ 100 Mq (ISmithetalL[2010cJ) . 

An amount of about 0.03-0.35 Mq of dust present in a circumstellar torus created by 
possible LBV-like mass-loss o f mass 3-35 Mq has be en proposed to be the origin of a mid- 
IR excess toward SN 2010jl ([Andrews et all 201 lah at ~ 90 d. From spectropolarimetry 
obtained at an earlier epoch (14 d), the p resence of a sign ificant dust mass along the line 
of sight to the progenitor was ruled out jPatat et al [201 ih , so the dust must lie in an in- 
clined torus. A warm Spitzer/IRAC survey of 68 Type Iln SNe detected between 1999 and 
2008 results in about 10 Type Iln SNe exhibiting late-time mid-IR emission caused by pre- 
existing dust hea ted through the interaction of the SN shock with the circumstellar medium 
jFox et a l[201lt). The progenitor mass-loss histories are consistent with those of LB Vs. 

[Smith et a] ([2003) measured the mass of a 19th century eruption from the well-known 
LBV 77 Car t o be about 12-20 M q. A dust mass of 0.4 ±0.1 Mq surrounding 77 Car was 
estimated bv [Gomezet"ail ( [2010l) who also estimated that > 40 Mq of gas has been ejected 
so far. SN 1961V was tentatively classified as an rj Car-like outburst with optically thick 
dust in a massive sh ell suggested to be present based on the fad ing of the light curve 
after around 4 years ( [Goodrich et aj^ I989[: Filippenko et"ai 19951). However, the nature 
of SN 1961V is contentious (e.g., [Stockdale et all [200l[ ; [Van Dvk et al [2002[ ; [Chu et all 
[2004[) . [Kochanek et al[ ( (201 0') suggest it to be a peculiar, but real SN which has experienced 
enhanced mass-loss prior to expl osion. Similar objects are SN 1954 J (.Smith et aL ,20011 ; 
[Van Dvk et aj [2005[) . SN 199 7bs ([Van Dvk et all \200( h SN 2000ch jWagner et all [2004h 
SN 2 002kg, and SN 2003gm ( [Weis and Bomansl [20051 ; [Maund et al [20061 ; [Van Dvk et all 
[20061) . 

The transients UGC 2773-OT and SN 2009ip ( [Smith et al [2010b[ ; [Folev et al [2010h 
were both LBV outbursts. The progenitor of SN 2009ip was serendipitously observed 10 yr 
prior to its outburst as an extremely luminous st ar and the mass was estimated to be about 
50-80 Mq ( [Smith et al[2010b[ ; [Folev et ai[2010l) . UGC 2773-OT was less luminous with a 
mass of > 25 Mq, but found in a ve ry dusty environment. Finally, a dusty nebula around the 
object HR Car ( [Umana et al[200^ consisting of amorphous silicates indicates that dust has 
formed during the LBV outburst. 

[Smith and Owocld[ ilOO &i deduced masses for the observed nebulae of several LB Vs and 
LBV candidates and concluded that an LBV giant eruption typically involves 10 Mq of ma- 
terial. The expansion veloci ties of such outbursts can be as high as 750 km s~^ as measured 
for 77 Car ( [Pavidsod . [l97l[) and up to 2000-3000 km s ' for SN 1961V ( [Goodrich et al 
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Table 3 Observed and derived properties of SNe 



SN 


SN Type 


Progenitor 


Mp [Mef 


tpe [d]'' 




Ti [K]^ 


Refs. 


2007od 


II-P 


SAGB 


-9.7-11 


300 


1.7 X lO-'' 


580 


1,2 










455 


1.9 X 10-" 


490 


2 










667 


1.8 X lO-'' 


600 


2 


2007it 


II-P 


— 


~ 16-27 


351 


1.6-7.3 X lO""* 


500 


3 










561 


7.0 X 10-5 


700 


3 










718 


8.0 X 10-5 


590 


3 










944 


4.6 X 10-5 


480 


3 


2006jc 


pec. Ibn 


LBV 


-40 


200 


6.9 X 10-5 


800 


4,5 










230 


3 X lO-* 


950 


6 


2005af 


II-P 


— 


— 


214 


4 X lO-* 


— 


4,7 


2004et 


II-P 


RSG 


9 


300 


3.9 X 10-5 


900 


4, 8,9 










464 


6.6 X 10-5 


650 


8 










795 


1.5 X 10-* 


450 


8 


2004dj 


II-P 


RSG 


12-20 


267-275 


0.3-2.0 X 10-5 


710, 186 


10-14 










500 


2.2 X 10-5 


650 


15 










652 


3.2 X 10-5 


610 


15 










859 


3.3 X 10-5 


570 


15 










849-883 


0.1-3.2 X 10-" 


530, 120 


10 










996 


5.0 X 10-5 


520 


15 










1006-1016 


0.1-7.6 X 10-" 


462, 110 


10 










1207 


> 1.0 X 10-" 


460 


15 










1236-1246 


0.1^.2 X 10-" 


424, 103 


10 










1393 


> 1.5 X 10-" 


430 


15 


2003gd 


II-P 


RSG 


-8 


499 


2.0-17 X 10-" 


480 


4, 16 










496 


4 X 10-5 


525 


17 










678 


2.7-20 X 10-^ 




16 


1999em 


II-P 


RSG 


15 


510 


- 10-" 


510 


4, 18 


1998S 


Iln 






360 


> 2 X 10-3 


1250 


4, 19 


1987A 


Il-pec 


BSG 


-20 


615 


3.7-31 X 10-5 


422 


4,20 










615 


2-13 X 10-" 
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775 


5.9-50 X 10-5 


307 


20 










775 


2-7.5 X 10-" 




21 










1144 


5 X 10-" 


150 


22 










6067 


1-20 X 10-" 


90-100 


23 










8467, 8564 


4-7 X 10-' 


17-23 


24 



References. (1) Insena et al (201 ll); (2) Andrews et al (2010); (3) Andrews et al (201 lb); (4) Smartt et aj 
(2009, and references therein); (5) ISakon et al (2009); (6) Mattila et al (2008a); (7) Kotak (2008); (8) 
iKotaketal (2009 ); (9) [Magu ire et al (2010); (10) Szalai et al (2011); (11) Mai'z-ADellaniz et al (2004); 
(12) Kotak etal (2005D; (13) Wangetal (2005); (14) Vinko et al (2009); (15) Meikle et al (2011); (16) 
Sugerman et al (2006); (17) Meikle et al (2007); (18) Elmhamdi et al (2003); (19) Pozzo et al (2004); (20) 
Wooden et al U993) ; (21) .Ercolano et al 12007.) ; (22) .Dwek et al ( .1992) ; (23) .Bouchet et al 12004) ; (24) 
Matsuura et ail <2011l) 

Notes. "Mp is the mass of the progenitor, ^fpe is the time past explosion. '^Mj is the inferred dust mass. '^T^ 
is the inferred dust temperature. 



Il989h . Dust formed in such LBV outbursts is likely to escape before the shock from the 
final SN explosion catches up with the dusty shell. 
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Fig. 2 Inferred amount of dust from SN and SNR observations at different (a) epochs and (h) temperatures 
(Tables|3]and|4). Filled circles represent observations of SNe at early and late epochs and open circles mark 
observations from SNRs with an age of several 100 yr The colours denote the temperature (Ta) of the dust 
and fpe is the time past explosion. The size of the symbols is scaled by the mass of the SN progenitor The 
horizontal dashed line represents an upper limit to the dust mass of 3 X 10"' M© at early epochs and is also 
consistent with the logarithmic average of the inferred dust masses in SNRs. 
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Table 4 Observed and derived properties of SNRs 



SNR 


SN Type 


Progenitor Mp [Mq]*"' 


t fvrl'' 




^ a l^^l 


Refs. 


Cas A 


lib 


WIS. ij — jyj 




1 1 X W '^ 


1 70 


1 2 








J ZD 


J.o X iU 




9 
z 










1 1 n-4 
— iU , — lU 


on 
JJU, vu 


a 
J 










^ V 1 n-^ 1 V T 

J X iU , i X iU 


zoo, zzo 


/I 










J X lU 


7Q S9 
/y, oZ 












< 1 5 












335 


2—5 4 X lO"*^ 


40—150 


5 








337 


- 1 


- 20 


7 








337 


6 X 10"^ 


- 35 


3 








337 


7 5 X lO""^ 


- 35 


9 


Kepler 


la /lb 


— ~8 


405 


1-2 X lO-"* 


107 


10, 11, 12<^ 








405 


5 X 10-"* 


75-95 


lo*- 








405 


0.1-1.2 


16, 88 


13 


B0540 


II-P 


— 15-25 


700-1100 


1-3 X 10-3 


50-65 


14-16 


Crab 


n-Por 


— 8-10 


950 


1-7 X 10-2 


45 


16-19 




ECSN 




950 


3-20 X 10-3 


50 


19 








952 


1-10 X 10-3 


74 


20 


1E0102 


Ib/Ic or 


— -30 


-1000 


1.4 X 10-2 


50-150 


21-23 




II-L/b 




-1000 


3 X 10-3 


70 


24 








-1000 


2 X 10-5 


145 


24 








-1000 


8 X 10-* 


120 


25 



Refere nces. (1) iKrause et all <2008l) ; (2) lArendt etaij < 19991) : f3 ) iDouvion et all <2001b l); (4) iHines et all 
20041) : (5) 'Wilson and Batrla' (2005): (6) Rho et al (2008): (7) Dunne et al (20091): (8) [Sibthorpe et a| 
20091): (9) Barlow etal (2010): (10) Blair et al (2007): (H) Reynolds et al (2001: (12)lDouvion et aj 
2001a): (13) Gomez et al (2009): (14) Reynolds ( 1985): (15) Williams et al (2008): (16) Chevalier (2006): 

(17) fNomoto et al (1982): (18) Kitaura et al (2006): (19) Green et al (200 4J: (20) Temim gt al (2006): 

(21) Blair etal (2000): (22) jChevaliej feOOSi) : (23) lRhoetai fe009l) : (24) ISandstrom et all (2008) : (25) 

Istanimirovic et al ( 200j) 

Notes. "Mp is the mass of the progenitor, ''fpe is the time past explosion. '^Mj is the infeired dust mass. '^T^ 
is the inferred dust temperature. ^The derived dust masses are attributed to circumstellar dust heated by the 
SN blast wave. 



5.3.3 Type la, lb, Ic and lib supernovae 



Significant amounts of dust from Ic or lb SNe has not been reported, and they are not cur- 
rently considered to be important sources of dust. 

A very clear non-detection of dust for a Ic SN was obtained by iHunter eta 1 ( l2009l) 
for SN 2007gr. Besides the peculiar Ibn SN 2006jc, the only proposed o ccurrence of dust 
formation for a lb SN is for SN 19901 at d ay ~ 250 fe lmhamdi et alLl2004) . The same seems 
to be the case for Type lib SNe. However. iKrause et all ( 120081) has identified the SN causing 
the SNR Cas A as a Type lib. Cas A is well studied in terms of dust (see Sect. 15.3.1b and 
represents the only example so far of a SN of this Type where dust has been reported. 

Cl avton et al ' ( 1997) disc uss the possibility of SiC grain format ion and growth in Type 
la SNe, but the latest models jNozawa et al, 201_13 and observations jBorkowski et al, 2006) 
of Type la SNe indicate that only little or no dust forms in the ejecta. Ilshihara et all ( 120 101) 
attributed the thermal dust emission in parts of the Type la SNR Tycho S N 1572 to a pos- 
sible S N shock interaction with ambient molecular clouds. On the contrary. |Tian and Leahvl 
l l20Ilh conclude from radio and X-ray observations that the remnant is isolated. The SNR 
Kepler possibly constitutes an exceptional case, with inferred FIR dust masses up to 1-3 Mq 
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jMorgan et alLl2003l : lGomez et a]Ll2009l) . However, MIR observations rather indicate a dust 
mass of a f ew 10"^ Mp of dust, sugg ested to arise f rom circumstellar dust heated by the SN 
blast wave jPouvion et all 2001a: Bl air et a i l2007h . Moreover, the classification of t he pro- 
genito r is debated. The first claim that Ke pler has its origin in a SN la were made bv lBaadj 
( Il943h which was later also supported bv lBlair et all ( l2007h . iBandieral lll987h suggested that 
the progenit or might have been a runaway star with strong winds. Further possibilities are 
discussed bv lRevnolds et all l l2007h . but a SN la event is favoured. Pertaining to the meagre 
evidence of dust from this Type of SN and its ambiguous nature and delay times, SNe la are 
likely not significant dust contributors in the early Universe. 



5.3.4 Summary of observational status 

The observational status of SNe and SNRs for which dust formation has been inferred is 
summarized in Tables[3]and|4l We have attempted to include all observed SNe and SNRs for 
which information about the derived dust mass is available. In addition, we present details 
about the Type of the observed SN, the nature and mass of the progenitor, the epoch of 
observation and the inferred dust temperature. Based on this we plot in Fig.[2]the observed 
dust yields from Tables [3] and |4] as a function of epoch (Fig.|2^) or temperature (Fig.|2]3). 
From Fig. [2^ it is evident that regardless of SN Type or progenitor mass, only hot dust at an 
amount below ~ 3 x 10"^ Mq is present at early epochs, i.e., less than about 2 x 10^ days 
past explosion At late epochs (later than ~ 5 x 10^ days past explosion) a large dispersion 
in the inferred dust masses is evident, spanning 7 orders of magnitude (from 10"^ Mq to 
1 Mq). One might speculate that the upper envelope indicates that, with aging of the SNe 
and SNRs, dust grains grow to larger sizes and the total amount of dust increases. However, 
the presence of only hot dust at early epochs might as well reflect an instrumental selection 
effect (see Sect. l5.2.7T l. At these epochs no observations at longer (FIR to (sub-)millimeter) 
wavelengths, sensitive to cold dust temperatures, have been accomplished so far. Thus, the 
presence of cold dust at early epochs cannot unambiguously be ruled out. 

Higher inferred dust masses appear to be related to cold dust. This is clear from Fig.|2j5 
which exhibits a conspicuous relation between the inferred dust mass and the temperature 
at which it is inferred. Independent of the epoch of observation or the progenitor mass, the 
amount of dust at lower temperatures is significantly larger than at warm to hot temperatures. 

Finally we point out that the progenitors of the observed SNe and SNRs, which are in 
the mass range of 8-30 Mq, eject a total mass of heavy elements relevant for dust formation 
of about 0.3-2 Mq in the SN explosion. Only stars more massive than ~ 15 M r eject an 
amount of heavy elements l arger th an 1 Mq (e.g., Wooslev and Weave; 1 19951 . hereafter 



IWW95h . lNo"moto et all ( l200 ^. herafter lN06l) . lEldridge et all ( l200 ^. herafter'ETOS')). Thus, the 
observationally inferred high dust masses in e.g., SN I987A (Table [3]l or Cas A (Table O 
likely necessitate very high dust formation efficiencies if of SN origin. The possible dust 
formation efficiencies and their uncertainties are discussed in following section. 



6 Dust production efficiency 

Based on the dust yields obtained from observations and theory, summarized in previous 
sections, we next discuss the efficiencies of massive stellar sources in producing dust from 
their available metals. 
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Fig. 3 Dust production efficiencies of AGB stars. Tlie efficiencies are based on dust yields from 
Iperrarotti and Gail (2006) and yields of heavy elements from van den Hoek and Groeneweaen ( 1997). The 
solid, dotted, dashed, dashed-dotted and dashed-dot-dotted curve are for metallicities of Z = 0.001, Z = 0.004, 
Z = 0.008, Z = 0.02, and Z = 0.04, respectively. The green line indicates the metallicity-averaged efficiency 
cagbCw) obtained as a straight average of the five black curves. The vertical dashed line marks the boundary 
of 3 Mq, below which AGB stars are not considered as dust contributors at high redshift. 



Following ' Gall et a l( l2011ah . the dust production efficiency e{m, z) per stellar mass and 
metallicity is defined as 

Mz{m, Z) 

where M^im, Z) is the mass of dust produced and released into the ISM, Mz(m, Z) is the total 
ejected mass of heavy elements relevant for dust condensation per star and m = A/»/Mo, 
where is the zero age main sequence mass. It is assumed that the amount of dust is the 
final mass, which has formed and possibly been processed through shock interactions. 



6. 1 Efficiencies of AGB stars 



The dust production efficiency, ep£}s(m,Z), for AGB stars in the mass range 3-7 Mq is cal- 
culated from theoretical values of Md (m,Z) and M7(m,Z ) . The amount of dust Mi(m, Z) 
is obtained from total dust yields of Ferrarotti and Gaill i 200q) . Stellar yields for AGB 
Stars have been calcula t ed by e.g.. lRenzini and Volil (Il981h.|] 
iKarakas and Lattanzid ( l2007h and recently by Karakad 1 20 



Marigol ( l200lhlHerwi"3 ( l2004h . 



However most of the mod- 



els do not provide yields covering the range of masses (3-8 Mq), elements or metallicities 
relevant for this investigation. For the sake of consistency, the amount of heavy el ements, 
Mz(m,z), is obtained from the the yields of Ivan den Hoek and GroenewegenI ( Il997t) cover- 
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ing a large grid of metallicities and stellar masses. The efficien cy eACjB(m,Z) is calculate d 
for four different metallicities in accordance with calculations bv lFerrarotti and Gaill ( l2006h . 

The results are presented in Fig. [3] It is evident that ej^QB(in,Z) decreases quite rapidly 
between 4 and 5 Mq, independently of the metallicity. AGB stars in the mass range 3-4 
Mq (i.e., C-stars) apparently are the most efficient dust producers. It can also be seen that 
at lower metallicities (Z < 0.008) these AGB stars are more efficient in condensing their 
available heavy elements into dust than at higher metallicities. The green thick curve in 
Fig. [3] illustrates the metallicity-averaged efficiency, eAGaim), for AGB stars. 



6.2 Efficiencies of CCSNe 

As highlighted in Sect. |5] there is a discrepancy between the derived SN dust yields from 
observations (resulting in low amounts of dust) and theory (predicting large dust masses). 
It is therefore of interest to determine plausible limits for the dust production efficiency of 
SNe based on the dust yields obtained from either approach. 

6.2.1 Maximum efficiency 

An upper limit to the SN dust producti on effi ciency can be ascertained using the mass and 
metallicity dependent dust yields from lTF()ll to deter mine th e mass of dust Mi(m,Z). The 
yields for the heavy elements Mzim, Z) are taken from lWW9i since these were also used by 
ItFOII . Th e effi ciency for Z = is derived from the dust yields of N03 and the total amount of 
metals of iNOfil - both yields are taken from unmixed grain models. The resulting efficiencies 
can be seen in Fig.|4l where we notice a clear decline of e(m, Z) with increasing progenitor 

mass. ^ ^ 

1 1';' efficiencies obtained bv lTFOll and InOBI at Z = differ significantly. For ItFOII , 
e{m, Z) decreases quite drastically for stars between 20-25 Mq, whereas it remains more flat 
in the models of N03. The maximum SN efficiency limit can be obtained by averaging the 
efficiencies obtained for each Z from these models over metallicity (the average efficiency 
for all stellar masses is obtained via rational spline interpolation and extrapolation into the 
mass regime of 8-12 Mq where no yields for heavy elements are available). This is sufficient 
to describe the observed tendencies and obtain an estimate of e(m). We will refer to this as 
the 'maximum' SN dust production efficiency emax('w), drawn as the dark blue curve with 
red crosses representing the averaged values in Fig.|4l 

6.2.2 High efficiency 

According to the predictions of lBianchi and Schneidej ( l2007h and lNozawa et all ( I2007l . l201oh . 
dust grain destruction takes place when a reverse shock penetrates the dust layer at timescales 
up to ~ IC* years past explosion. This leads to a significant (up to 10 0%) reduction of the 
dust formed, depending on the ISM density and grain size. For example. lNozawa et ail( l2007h 
have shown that large gr ains in contrast to small grains remain relatively unaffected by the 
reverse shock. Following lBianchi and Schneided ( 20071) , the possibility of grain destruction 
can be accounted for by applying a reduction of 93% to ej^nxim). The resulting reduced 
efficiency still represents a rather high dust production efficiency in comparison to what is 
derived from SN observations. Hence, this will be referred to as the 'high' SN dust efficiency 
ehigh(w?) = 0.07en,ax('w)- 
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Fig. 4 Dust pr oduction efficiencies of mas sive stars. Upper curves: Efficiencies calculated from SN dust 
yields of iTFOll and metal yields of WW95. The solid (thin), dotted, dashed and dashed-dotted curve are 
for metallicities of Z = 0, Z = 0.0001, Z = 0.01 and Z = 0.02, respectively. The dashed-dot-dotted curve 
represents the efficiency at Z = derived from dust yields of N03 and metal yields from N06. The blue thick 
curve represents the averaged 'maximum' SN efficiency fniaxC"') and the black thick solid curve represents the 
fitted efficiency £^.^^{m). The thick violet curve represents the 'high' SN efficiency ehishC'")- Lower symbols: 
Efficiencies derived from the averaged observed dust amount of 3 X 10"^^ Mq of SNRs and the SN metal 
yields of WW95 (stars), NOfi (crosses) and ET08 (triangles) for solar metallicity. The thick cyan curve is the 
averaged 'low' SN efficiency eiawitn) and the black thick solid curve the fitted efficiency e^^^{m). The left 
green curve represents the averaged AGB efficiency eAGB('") (see also Fig. [3). The vertical fines mark the 
range of AGB stars between 3-8 Mq and SNe more massive than 8 Mq. 
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We note that either CmaxC™) or ehigh(w') might also be interpreted as the result of longer 
timescale dust grain growth (see Fig. |2]l in the SNR itself. The 'maximum' emax(w') pre- 
supposes that dust destruction through shock interactions is inefficient. The 'high' SN ef- 
ficiency, ehish('n), could also be the result of smaller or no destruction, depending on how 
much dust would initially have formed before a possible shock interaction. 



6.2.3 Low efficiency 

The lowest feasible limit for the efficiency of SNe dust production is generated based on 
observed dust yields from the SNRs Cas A, B0540-69.3, Crab nebula, and 1E0102.2-7219 
at temperatures between 50-100 K (see Table|4]i. The inferred amount of dust Ma(in,Z) is 
taken to be 3 x 10"^ Mq and is applied to SNe in the mass interval 8^0 Mq^ 

For the mass of heavy elements Mz(m, Z) the yields of lWW95l InO^ and lETOSi are used. 
The metallicity of most SN progenitors given in Tables [3] and |4] is estimated to be between 
around solar (Z = 0.02) or LMC-like (Z = 0.008) (Smarttetal, 2009). Solar metaUicity 
for all SNe in the mass range of 8^0 Mq can therefore be assumed. The metal yields 
Mz{m, Z) are also evaluated for Z = Zp. To obta i n the low SN efficiency limit we average the 
efficiencies obtained using the yields of lWW95l. lN0^ and ET08. The same interpolation and 
extrapolation scheme as for e,„ax(m) is applied and the resulting average efficiency appears 
as the cyan curve with average values indicated as red crosses in Fig.|4] 

The resulting averaged dust production efficiency only depends on the stellar mass. We 
will refer to this as the 'low' SN efficiency eiow('«)- Interestingly, also eiow('w) features a 
declining tendency with increasing stellar mass, similar to emax('w). 

There are two possible interpretations of this limit. The amount of dust produced by 
SNe could be similar to the low observed amount of dust at early epochs and this rather 
low amount of dust does n ot significantly grow on longer timescales. This might be the case 
for the SNR B0540-69.3 dWilHams eta]li2008h . Alternatively, eio„(m) may be the result of 
potential dust destruction of larger amounts of dust from shock interactions. 



6.3 Analytical approximations 

To illustrate the general trends of different e{m) we provide simple analytical fits to the 
derived averaged efficiencies of AGB stars and SNe. 

One notices from Fig. |4] that there might be a smooth connection of the efficiencies, 
eiov/(m), between high-mass AGB stars and low-mass SNe. An adequate approximation cov- 
ering all stars between 3^0 Mq is a power law for e(m), 

eiow(™) = « "''^ + 3 < m < 40, (3) 

with a = 15, /? = 3.25, and c = 2.8 x lO""*. The negative slope reflects the decreasing ef- 
ficiency of stars with increasing mass to release the produced dust grains into the ISM. It 
also illustrates that AGB stars in this case are more efficient, closely followed by the low- 
mass SNe. While e^^^(m) drops by roughly three orders of magnitude in the 3^0 Mq mass 
range, the rather steep decline for stars between ~ 3-12 Mq over approximately two orders 
of magnitude is noteworthy. We also note that although e^^^{m) provides a fairly good ap- 
proximation to eiow(wj), it does not capture the strong preference for 3^ Mq stars over 5-7 
Mq stars (see Fig.[3]l. 
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Fig. 5 The five IMFs considered. The solid, dotted, dashed, dashed-dotted, and dashed-dot-dotted curves 
represent the Salpeter, mass-heavy, top-heavy, Larson 1 and Larson 2 IMF, respectively. The dark grey area 
signifies the critical mass range of 8-12 M© and the light grey area signifies SNe between 12-20 M©. The 
vertical dashed lines mark the limits of the range in stellar mass considered (between 3 M© and 40 M©). 



The 'maximum' SN dust formation efficiencies are better approximated by an exponen- 
tial function. 



-(m/jiio) 



8 < m < 40, 



(4) 



with a = 1.2 and mo = 13. Comparing the efficiency eAGs(m) of AGB stars to e^n^im), we 
find no possibility for a smooth connection. In this case, stars between 8-12 Mq are the most 
efficient dust producers. The general decline of e^ax("^) for stars between 8^0 Mq is about 
an order of magnitude, comparable to the drop of eAGB('w) from a 4 Mq to a 6 Mq AGB star. 
The resulting fits are shown in Fig.|4]as black solid curves. 



7 Stellar dust productivity 

Having reviewed the dust production efficiency of single massive stars, we next discuss the 
dust productivity of these massive stars in a galaxy. Besides the dust production efficiency, 
the total amount of dust produced in a galaxy depends on the star-formation rate (SFR), i//(t) 
and the IMF, 4>{m). 



7.1 The initial mass function 



The IMF is an important parameter influencing the evolution of dust, gas and metals in a 
galaxy. It determines the mass distribution of a population of stars with a certain ZAMS 
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mass. The IMF was first proposed bv lSalpeteil ( ll955h . derived for Galactic field stars. Orig- 
inally, the IMF was not a power law but composed of a logarithmic slope of about -1.7 
for stars below 1 Mq and -1.2 for stars between 1-10 Mq. It was suggested that a power 
law with slope -1.35 applied to the entire mass range is appropriate, but strictly speaking 
it is only valid for stars between 0.4-10 Mq. Nevertheless, the Salpeter IMF is still often 
applied to more extended m ass ranges (e.g., 0.1-100 M©) (for detailed reviews, see, e.g., 
IScalol l2005l : IChabriej , l2005h . Later studies have shown that the IMF flattens for stars be- 
low 0.5 Mq and significantly declines in the mass regime < 0.1 Mq (e.g., Kroup a, 2002.; 



IChabrieiLl2003allbl) A steeper decline for intermediate mass stars has also been suggested 
llScald.ll986l . 119981). A characteristic mass has been defined such that half the initial mass 



goes into stars with masses lower than the characteristic mass and half into stars more mas- 
sive. The characteristic mass describes the mass at which stars are preferentially formed . 
For th e field star IMFs described above the characteristic mass is about 1 Mq (e.g.. lLarsonl . 
l2006h . 

A fundamental debate regarding the IMF is whether there is a systematic variation of 
the IMF with som e physical conditions of star formation or whether it is universal (see 
iBastian et ail2O10L for a review). Several theoretical approaches suggest non-universality. 
Systematic changes of the IMF leading to a shift of the characteristic mass towards higher 
stella r masse s are fo und in star-formin g environments with increased ambient temper atures 
(e.g.. lLarsorj.ll998h or high-densities jMurrav and Linl. Il996l : iKrumholz etaj l2010h . The 



absence of a variation of the characteristic mass due to a change of the equation of state as 
a result of dust processes is suggested by Bonnell e t al (2007). 

Usually any IMF with a characteristic mass shifted towards high stellar masses re- 
sulting in an overabundance of high-mass stars is referred to as 'top-heavy' IMF. The 
possibility of a top-heavy IMF in low- metallicity envirormients and in pa rticular in the 
early Universe was suggested already by Schwarzschild and Spitzer ('1953') and plausible 
evidence is extensively discussed in, e.g.. iLarson tl998i) and .Tumlinsoa (,2006) . Further- 
more, indirect and direct evidence for a t op-heavy IMF has been found in various sys- 
tems such as e.g., s tarburst galaxies (ee.g., Rieke et al, 1993; Doane and MathewsL [l993l ; 
iDabringhauseii et al 2009i) , disturbed galaxies (.Habergham et al, 2010), and sub-mill imeter 



galaxies (e.g., iBaugh etaj l2005l ; iNagashima et aj I2OO5I ; iMichalowski et all 1201 Obi) . Evi 



dence for IMF variations towards higher stellar masses also comes from observations of 
the Galactic Center region and Galact i c globular cluster s (e.g., D'Antona and Caloi, 20041; 



Ballero et a i l2007l ; lManess et a]ll2007l ; lBartko et alll2oTd) and star clusters (e.g.. lSmith and Gallaghej , 



20011) 



However, any reference to a top-heavy IMF must be assessed critically. Usually the 
degree to which the IMF is 'top-heavy' varies significantly among different studies. Differ- 
ences can be due to extreme assumptions of the exponent for the power law IMFs (i.e., a = 
0) or due to different characteristic masses in log-normal IMFs, but may also be due to the 
assumed mass interval of the IMF. 

In models of dust e volution in galaxies and high-z QSOs (e.g.. [Morgan and Edmundj , 
l2003l ; lDweke"tIill2007l) a Salpeter IMF is often used. Considering the evidence for an IMF 
different from the commonly adopted Salpeter IMF, it is of interest to investigate the impli- 
cation on the dust productivity of five different IMFs (Table[5]l. 

The IMF is normalized in the mass interval [mi, m2] such that 



m (p{m) dm = 1, (5) 
where mi and m2 are the lower and upper limits of the IMF given in Table |5] 
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Table 5 IMF parameters 



IMF 


a 


nil 


1112 


"Jch 


Salpeter 


1.35 


0.1 


100 




Mass-heavy 


1.35 


1.0 


100 




Top-heavy 


0.5 


0.1 


100 




Larson 1 


1.35 


0.1 


100 


0.35 


Larson 2 


1.35 


0.1 


100 


10.0 



The power law IMFs (Salpeter, mass heavy a nd top heavy) have the form (p(m) oc ' 
while the log-normal Larson IMFs tLarsonLri99&) are given as (f)(m) oc m"^""^'' exp(-mch/'«), 
where nich is the characteristic mass. The 'top-heavy' IMF is characterized by a flatter slope 
than the Salpeter IMF. The 'mass-heavy' IMF has a similar slope as the Salpeter IMF but 
the formation of stars with stellar masses below 1 Mq is suppressed leading to the formation 
of more stars in the mass interval [mi, WI2] compared to the Salpeter IMF. The log-normal 
IMFs have the same slope as the Salpeter IMF in the high mass tail of the IMF, but flatten 
or decline for masses below the characteristic mass. The 'Larson 1' is closest to a Salpeter 
IMF while the 'Larson 2' IMF is biased towards higher stellar masses and can be referred 
to as a 'top-heavy' IMF. 

In Fig. [5] we plot the IMFs considered. From the shape of the curves for (/)(m) it is 
evident that the majority of the stars relevant for dust formation are formed in the mass 
range of 3-8 Mq for all IMFs, followed by stars between 8-12 Mq. Note that this includes 
the critical mass range of 8-10 Mq (see Sect. 12.41 ). For SNe between 10-12 Mq dust yields 
or metal yields are uncertain or unavailable, leading to uncertainties in the dust production 
efficiency. For the purposes of the following discussion, we therefore extend the previously 
defined critical mass range up to 12 Mq. 



7.2 Dust productivity 

The amount of dust produced per star, Ma{m), is calculated from the dust formation ef- 
ficiencies as Mii{m) = Mz{m) e(m). The yields of heavy elem ents Mz(m) are taken from 
IWW95l (for SNe) and Ivan den Hoek and Groenewegenl ( Il997h (for AGBs). We study two 
cases, Z = Zq and Z = O.OIZ©. 

To quantify the effect of the various IMFs and the dust production efficiencies on the 
total dust contribution from AGBs and SNe we define the total dust productivity of all stars 
in the mass interval [m^, mj/] as 

r'"" Mz{m) 
fiD = <p{m) — — e(m) dm. (6) 

JniL "^O 

The lower and upper mass limits, mi and my, deliniate the interval 3^0 Mq, which will be 
further divided into the AGE star range 3-8 Mq, and the SN ranges 8-12 Mq, 12-20 Mq 
and 20-40 Mq. The total dust productivity, /io, depends on the IMF, the efficiency and the 
metal yields through the integrand fd('«) = 0(m) (Mz(»i)/Mo) e(m), which is the specific 
dust productivity. 

The calculated amount of dust M^im) for each efficiency case (see Sect. l6.2l i is presented 
in Fig. [6] We note that the amount of dust produced per AGB star is between the values 
of Mi{m) for SNe with e\o„{m) and ehish('«). Regarding the quality of the fitting functions 
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Fig. 6 Dust yields for AGB stai's and SNe calculated for e,„ax('JO (dark blue curve), ehiah('") (violet curve), 
eio„(m) (cyan curve), e*^(m) and e^^^(m) (black curves) as well as for AGB stars (green curve). Filled circles 
represent observed dust yields for different SNe at different temperatures. The dark grey zone corresponds to 
the critical mass range (8-12 M©) and the light grey region con'esponds to the approximate mass range for 
Type II-P SNe. 



(Eq. |4ll, using e^^^(m) for all stars in the range 3^0 Mq results in lower dust yields for 
stars between 6-7 Mq and a significant overestimate of for stars in the critical mass 

range 8-10 Mq, relative to using £iQ^(m). Using £^^{m) for the 'maximum' SN efficiency is 
consistent with using emaxini) 

For comparison we plot the highest inferred dust yields from the observed SNe listed in 
Tables[3]and|4] The two upper values for Cas A jRho et a]ll200i : lDunne et a]ll2009h and the 
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Fig. 7 Specific dust productivity, (a) Specific dust productivity ^d(m) of stellar masses calculated for 
(dark blue curve), ehigh(m) (violet curve), eiow('«) (cyan curve), ejj,ax('") and €^^^{m) (black curves) as well as 
for AGB stars (green curve). The dark grey zone corresponds to the critical mass range (8-12 Mq) and the 
light grey region corresponds to the approximate mass range for Type II-P SNe. The solid, dotted, dashed, 
dashed-dotted, and dashed-dot-dotted curves represent the Salpeter, mass-heavy, top-heavy, Larson 1 and 
Larson 2 IMFs, respectively. 
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upper value of SN 1987 A jMatsuura et all l20Tlh at low temperature match the dust yields 
calculated using ema\{ni) or ehigh(ffj). Dust masses for SNe calculated using eiow(«i) are also in 
good agreement with the observed dust yields from several SNRs. We also plot the highest 
observationally derived dust yield for the Kepler remnant (see Table |4] and discussion in 
Sect. l5.3.3] |. The inferred dust yields from the Kepler remnant may not be representative of 
dust from CCSNe (the SN has been suggested to be a Type la). However, in view of the 
general uncertainty about dust from stars with such a progenitor mass, this SNR provides an 
interesting benchmark. 

Fig.[7]shows the specific dust productivity fd('w) for all e(m) and the various considered 
IMFs. The slopes of ^d(m) exhibit a declining trend with increasing stellar mass regardless 
of the choice of e(m). SNe with masses between 30^0 Mq are ~ 10 times less productive 
than SNe with masses between 8-12 Mq. For AGB stars ^d{ni) decreases steeply between 
3-7 Mq, resulting in about an order of magnitude lower value for the higher mass AGB stars. 
The most productive AGB stars therefore are 3-4 Mq stars, partly reflecting their higher dust 
production efficiencies eAGB('w) (see Figs. |4] and IS}. 

A Larson 2 IMF exhibits the highest dust productivity for SNe between 8^0 Mq, inde- 
pendently of the dust production efficiency, while the lowest specific productivity is obtained 
for a Salpeter IMF. The difference in fd('«) between either a Larson 2 or a top-heavy IMF 
and the Salpeter IMF is larger for the more massive stars (~ 30-40 Mq). For AGB stars, the 
largest sensitivity to the IMF occurs for 3^ Mq stars which exhibits the largest difference 
in ^d(m) for a mass-heavy IMF (highest value) vs. a Larson 2 IMF (lowest value). 

The total dust productivity jjo of AGB stars and SNe, subdivided into 3 mass ranges, 
is presented in Fig.[8l For a 'low' SN efficiency, eiow(»j), the total amount of dust produced 
is almost exclusively manufactured by AGB stars. Dust production by SNe in this case is 
negligible for all considered IMFs. 

The total dust productivity is increased as soon as SNe are assumed to produce dust with 
the 'high' SN efficiency ehigh(m). For a Salpeter, Larson 1 and mass-heavy IMF, AGB stars 
still dominate the dust production whereas for a top-heavy or Larson 2 IMF, SNe are the 
prime dust producers. 

In case of the 'maximum' SN efficiency, emnx(m), dust is primarily manufactured by 
SNe and the dust supply from AGB stars is negligible. For this efficiency the amount of dust 
produced by SNe is roughly 5-10 times higher than for eiow('w) and ehigh(m), depending on 
the IMF. We find that for the IMFs favouring lower mass stars, the three SN mass ranges 
(8-12 Mq, 12-20 Mq, 20^0 Mq) are nearly equally important. While there is considerable 
uncertainty about dust production channels from stars between 30^0 Mq, the analysis in- 
dicates that this mass range is the least significant range. Hence, SN dust production is in 
general dominated by stars in the mass range 8-20 Mq with an almost equal contribution 
from stars between 8-12 Mq and 12-20 Mq. 

For the calculations with yields for heavy elements at a metallicity of Z = 0.01 Zq we 
find the same tendencies. This indicates that these relations most likely also apply to high-z 
galaxies. 

7.3 An example 

This simple formalism allows us to address the origin of large dust masses in QSOs at high 
redshift, as discussed in more detail in Sect. 18.11 In the following we estimate whether the 
derived dust productiv ities may be sufficient to account for the 2-7 x 10^ Mg of du st inferred 
in QSOs at z > 6 (e.g. lBertoldi et alll2003l ; lRobson et all[2004lBeelen et alll2006l) . 
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Fig. 8 Total dust productivity of AGB stars and SNe for different IMFs and SN dust production efficiencies. 
Tlie heiglit of tlie bars represents tlie total dust productivity of stars in the mass range 3^0 M©. The contri- 
bution from AGB stars is marked in green, the SN mass ranges of 8-12 Mq, 12-20 Mq and 20-40 are 
the dark grey shaded, Hght grey shaded and solid blue areas. The contribution from stars in the mass range of 
30—40 Mq is the area from the bottom to the black solid line in the bar. The letters S, LI, M, T, L2 stand for 
the Salpeter, Larson 1, mass-heavy, top-heavy and Larson 2 IMFs, respectively. The red dashed line marks 



the minimum estimated dust productivity of /jd ■ 
z>6. 



10 required to account for high dust masses in QSOs at 



We assume a minimum required dust mass of Mq = 2x10^ Mq and a maximum available 
time span of At = 400 Myr for building up this amount of dust. The minimum required 
average amount of dust produced per unit time is expressed as the dust production rate in 
this period, Ro = Mo/At = 0.5 Mq yr~' = //d iA(0- We assume a high constant average SFR 
i//{t) = 500 Mq yr~' , based on derived SFRs from ob served hi gh-z QSOs ranging fr om 100- 
3000 Mq yr"' (e.g. lBertoldi et alll2003l : lDwek et all 2007; Ri echers et al 2009; W ang et all 
bOldi . With these assumptions, all cases for which yUo ^ 10"^ can be excluded (see Fig.[8]l. 

For the 'low' SN dust production efficiency, e\ow(m), none of the IMFs gives a suffi- 
ciently high dust productivity. Only a mass-heavy IMF is close to the limit. Moreover, the 
long lifetimes of 3^ Mq AGB stars (see Fig. O, which dominate the AGB dust production 
(see Sect. |7] and Fig.[6ll, is problematic. These stars will start contributing with a delay of 
more than ~ 200 Myr so will produce dust for only approximately half the time of the as- 
sumed maximum time span of 400 Myr. Thus, a 'low' SN dust production efficiency appears 
to be insufficient to account for the dust at high redshift. 

In case of a 'high' SN efficiency, ehigh(m), the majority of the IMFs might lead to a suf- 
ficiently high dust productivity. Due to their short lifetimes, SNe can be assumed to release 
dust immediately after formation. Thus, SNe dominate the dust production for a Larson 2 
or top heavy IMF. Taking into account the reduction of the AGB dust contribution due to 
the long lifetimes of these stars, the Salpeter or Larson 1 IMFs most likely do not lead to 
sufficiently large amounts of dust at high-z. 

For a 'maximum' SN efficiency en,ax('«) the total dust production rates Ro of 3-18 Mq 
yr"' are achieved primarily through SN dust production. This leads to possible dust masses 
in excess of 10' M© produced in high-z systems, even for significantly lower star formation 
rates than assumed in our scenario. 

We note that it is unclear if a high SFR can be sustained over 400 Myr. In fact, the 
very high derived SFRs (> 1000 Mq yr~') are attributed to shorter (< 10^ yr ) durations of 
the starburst (e.g. lBertoldi et alll2003l ; lDwek et alll2007l ; lRiechers et all 12003) . Assuming a 
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SFR ip(t) = 1000 Mo yr"' and a At = 200 Myr leads to the same dust productivity /lo = 10"^ 
as discussed above, although the AGB star contribution would be even more suppressed due 
to the long lifetimes relative to At. 



8 Dust at high redshift 



From SC UBA, MAMBO, MAMBO-2 and VLA surveys of brigh t high-z Q SOs at 4 < z ^ 



6.4 (e.g..lCarim et alLl200Ial;IOmont 



003; I 



200lUlsaak et alll2002l;rB ertoldi and Coxl [20o3 : 
Beelen et a 11200^ very high 
Mq yr~' have 



iBertoldi et al l2003l : [Priddev et all 120031 : iRobson et all 12004 
dust masses of more than 10* Mq and star formation rates of more than 10^ 
been inferred from the measured sub-millimeter fluxes. 

Theoretically, it has been proven difficult to explain the origin of these dust masses in 
QSOs at z > 6, despite some attempts (e.g.. Dwek et al. 2007; Vah ante e t al, 2009 ; Pipi no et a^ 
|201 lllDwek and Cherchnefd.l201 iLlGail et ail201 l4 d lMattssoil(20rii ; l\^liante et ailbol ih . 



At this redshift the timescale available to build up large dust masses is short, which limits 
the possible options for sources of dust (see Sect.[2l Fig. [Hi. 

As a consequence, massive stars have been strongly favoured, although the actual dust 
production by massive stars (see above discussions in Sect. |4] and IS} is afflicted with large 
uncertainties and other sources may play an important role as well. 

Nevertheless, features in the extinction curves of various objects at high redshift have 
been attributed to dust of SN origin. For example, the extinction curve inferred for the QSO 
SDSS J1048+46 at z = 6.2 possesses a characteristic plateau at aro und 1700-3000 A, in- 
terpreted as arising from amorphous carbon and magnetite SN dust (lMaiolinoetalll2004) 
based on the Itfoi I models. A similar feature has been reported for the afterglow of GRB 
071025 at z ~ 5 jPerlev eTal \20l(h . Also less conspicuous features in extinction curves 
(notably flatter UV slopes than the SMC extinction curve) have been interpreted as evidence 
for dust from SNe. Several QSOs (e.g., Gallerani et al, 2010) tu rned o ut to be best fitted with 
a contribution from extinction curves for SN-like dust ( iHirashita et alll2008h . The young in- 
frared galaxy SST J 1 604+4304 at z ~ 1 ha s also been proposed to be best fitted with a SN 
extinction curve from lHirashita et all ( l2008h ( iKawara et allbOllh . 

The observationally derived dust masses and SFRs for high-z galaxies are naturally un- 
certain. Below we therefore briefly review the basic concepts and caveats in deriving dust 
masses and SFRs from observations. Next we summarize the theoretical models aiming to 
explain the observations at high redshift. 



8.1 Inferring physical properties of high-z QSOs 
8.1.1 Dust mass 

Based on the method discussed bv lHildebranldigSSh . the dust mass is determined from the 
sub-millimeter flux density observed at frequency Vg = vj{l -H z) as 

Mi = —-^ , (7) 

where is the rest-frame frequency (see also Eq. [T]l. For definition of the parameters we 
refer to Sect.lO 
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While there are uncertainties related to the cosmology entering the luminosity distance, 
Z)l> the main uncertainties in deriving the dust mass from observations are given by and 
Kdivr), vvhich is usually parametrized as 



KiiVr) = Kiivo) — , (8) 

where /? is the emissivity index. From the above formalism it is clear that A-ti(v) significantly 
depends on dust properties such as /?, the grain radius a and the grain density p (because 
= (3/4)2(v)/(ap)). None of these properties are well known. 
The dust absorption coefficient at the critical frequency (wavelength), vn = 2.4 THz 
(A{) =125 yum), at which a source becomes optically thin, was determined by Hildebrand' 
1 11983,) to be Kg ivo) = 18.75 cm^ g ~'. For values of the absorp tion coefficient oth er than that 
determined bv lHildebran3lll983l) . we refer to a summary of lAlton et all t004. their Table 



4). While Kii(v) increases from FIR to sub-millimeter wavelengths l DraineLll990h . it should 



be noted that even for similar wavelengths the inferred values for K^iv) often vary by an 
order of magnitude. 

Another ambiguous parameter is the emissivity index It has been found that is 
dependent on the wavelength (or frequency) and increases with increasing wavelength. For 
A < 200 yum the emissivity index ~ I and for A > 1000 yum, jS ~ 2 (e.g., lErickson et all 
1 1 98 iL ISchwartzL fl 982!) . However, fi might also depend on the dust comp osition, the grain 
size an d possibly also the temperature. For a detailed discussion we refer to lPunne and EalesI 
( |200I[) and references therein. The emissivity index /? as well as the dust tem perature, Tj, 
can be determined by fitting the spectral energy distribution (SED). According to lHildebrand 
lll983h the flux density, S (v), is defined as 

S{v) = QiQ{v)B{v,Ti\ (9) 

where Q^, = N{o-a/Dj) is the solid angle subtended by the dust source in the sky, with the 
number of spherical grains, each of cross section o-^. For high-z objects the SEDs are fitted 
in the rest-frame and Eq.|9]needs to be modified accordingly. 

For a simultaneous determination of Td and /3 many flux measurements at different wave- 
lengths are necessary. This however is often not possible f or hig h-z objects and values for 
either T^i or /? are simply assumed. iPriddev and McMahonI ( l200ll) found that the composite 
SED of a sample of QSOs at z > 4 are best fitted w i th a si ngle temperature of ~ 40 K 
and an emissivity index /J ~ 1.95, while lHughes et all l ll997h and lBenfordetai ' ( 1999) found 
Td ~ 50 K and /J ~ 1 .5 for high-z objects. From a study similar to IPriddev and M cMahonI 
l l200lh . but with a larger sample of high-z QSOs (1.8 < z < 6,4), Beelen et al (200^ obtain 
a higher temperature ~ 47 K but a lower ~ 1.6 for a combined SED of all QSOs (see 
Fig.©. 

The SED can in principle be fitted using either a single temperature m odel (as described 
above) or a two-temperature component model as accomplished by e.g.. iDunne and EalesI 
(l200ll) .rviahakis et al (2005) or Ivison et a], C20I0i) . For a two-component model the equa- 
tion for the dust mass can be expressed as 



5(yo)/)? 

Md 



(1 +z)Ki{v,) 



B{v,, T-w) B(v,, T,) 



(10) 



where 7V„ and represent the mass fractions of the warm and cold components. While 
the uncertainties in deriving the dust mass from observations are normally large, it has 
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been found that using a two-component dust model, the derived dust masses are usually 
a factor of ~ 2 higher than what can be ob tained from a single temperature model (e.g., 
iDuime and Ealesll200lhrviahakis et allbOOSh due to the larger amount of cold dust. 




Fig. 9 Combined SED, in the rest frame, of high-z QSOs taken from 'Beelen et al' I'lOOfT). Tlie plot co mprise s 
data from sources discussed in Benford et al (J^ggJJiilriddey and McMahon 1 2001) and Beelen et"ail i200d) . 
The mean FIR data points are best fitted with a graybody of temperature Td ~ 47 K and emissivity index 
/3 ~ 1.6. 



8.1.2 Star-formation rate 



The SFR of a galaxy, ijj, can be related to its dust continuum spectrum through the FIR 
luminosity, LpiR, as 

iA = 5MF5sB(£-FiR/10'"Lo)Moyr-i (11) 

(e.g., iGallagher et all [l98i; iThronson and TelescoL [losi; lOmont et all l200lh . Here 5mf = 
(^?FiR/Myr)"'(M/M0)/(L/L0), where zlfpiR accounts for the (assumed) duration of the star- 
burst, and MIL is the mass-to-luminosity ratio, which is determined from an assumed IMF. 
5sB is the fraction of the FIR emission due to dust heated by the starburst. 

The FIR luminosity can be obtained (see e.g.. lYun and Carillill2002h by integrating the 
SED (Eq.|9]l over the emitting area Q^. and the corresponding frequency range. 



JQi J 



(12) 
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Alternatively, once the dust mass is known, LpiR can be obtained by integrating the SED 
using Eq. (|7]l, 

Lfir = 47rMd J Ka(v)B(v, n)dv, (13) 
which emphasizes the relation between the FIR luminosity and the dust mass. For bright 



high- z objects, FIR luminosities of the ord er of 10'^ ' " ^ Lp are usual l y derived (e.g., Omont etalL 
200ll:[Berto ldi and Cox, 2002; De Bre uck et alll2003l:lRobson et alll2004lBeelen et alll2006l : 



Wang et alil20ldllvison et alil20iq : lLeipski et alll201ofy 



Evidently, the calculated value of the SFR sensitively depends on the as sumed IMF 
throug h ^ Mp. In most cases a Salpeter I MF is assumed, but, as pointed out by IPwek et all 
( l2007h and lDwek and CherchnefdteOUh . the IMF constitutes one of the major uncertainties. 
For e xample, the derived SFR of ~ 3400 Mq yr"' for QSO SDSS J1148+5251 CFan etall 
l2003h . using a Salpeter IMF, decreases to about 380 Mq yr ' for a top-heavy IMF. The range 
of the SFR in some high-z objects might therefore be between 10^"'* Mq yr"^ 

Another critical parameter is the ass umption of the durati o n of the starburst, AfyiR. Com- 
monl y either values of 6mf ~ 0.8-2.1 dScoville and Youn3 . 1 19831 : Ixhronson and Telescd . 
Il986h or simply Smf = 1 are adopted. However, these values have b een derived using a 
Salpe ter IMF and an assu med starburst age of for example /IfpiR = 2 Myr jXhronson and Telescol 
Il986l) . As pointed out bv lOmont et all ( 1200 ih , these assumptions might in fact be inappro- 
priate for massive starbursts in high-z galaxies. Consider ing a continuous st arburst of 100 
Myr and a Salpeter IMF with different low mass cutoffs, lOmont et all ( l200lh derive ^mf ~ 
1.2-3.8. Assuming a flat IMF ( a = 1) at lo w masses and zl/piR = 10-100 Myr results in Smf 
~ 0.8-2, similar to the values of iThronson and Telesco (1986). 

Regarding the fraction of the FIR emission heated by the starburst, 6sb, it is unknown 
whether the FIR luminosity arises solely from the starburst and if the entire stellar radiation 
is absorbed and re-emitted by warm dust or whether heating by an active galactic nucleus 
(AGN) must to be taken into account. The common view is that the heating source is the 
starburst and a contribution of the AGN is usually neglected, thus ^sb is set to 1. Taking a 
contribution of the AGN into account would result in a smaller a mount of dust and l ower 
star formation rat es. For a more detailed discussion we refer to, e.g., IOmontetall ( l200ll) and 
Ilsaaketall ( l2003) . 



8.2 Theoretical models 

To address the issue of the inferred large dust masses in high-z galaxies, chemical evolution 
models so far have been the preferred approach for following the temporal progression of 
the physical properties of a galaxy. The range of applications is large, e.g., the models can 
be used to investigate the temporal evolution of the abundance of different elements and 
dust influenced by formation and destruction processes, the abundance distribution of ele- 
ments, stellar masses, the metallicity, SFR and other physical properties. The models are 
mainly regulated by the interplay between processes such as star formation, gas and dust 
flows, stellar feedback, and the considered dust destruction and growth processes. Another 
important input is the IMF. However, most of the processes governing the models are uncer- 
tain and various simplifications must be made. For a profound review on chemical evolution 
models see Tinslev ( 1980), Dwek ( 1998), Dwek et al (2009) or Piovan et al (201 la b). Over 
the past years several dust evolution models addressing the above presented issue have been 
developed. Below we discuss the main findings. 
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The first atte mpt to explain the l arge derived dust m asses in the z = 6.4 QSO Jl 148+525 1 
(lFanetalLl2003h was carried out bv iDwek et a who concluded that at least 1 of 

dust per SN is required to account for the observed dust mass in this QSO (if SNe are the 
only sources of dust). The total mass of the QSO host galaxy is co nsidered to be abou t 5 x 
10^ ° Me to match the su ggested dynamical mass for this QSO (e.g. lWalter etail2004h . 

IVaUante et all ll2009l) include AGB stars in their models and claim that 10* of dust 
can predominantly be produced by AGB stars in QSO Jl 148+5251. The considered mass of 
the galaxy of about 1 x 10^^ Mq in the study, however, exceeds the plausible dynamical mass 
derived from observations by more than about an order of magnitude. The model includes 
a star-formation history resulting from a hierarchical galaxy merger tree scenario ihi et all 
I2OO7) and neglects gas in- and outflows. Star formation commenced at z = 15, resulting in 
abo ut 550 Myr availab le for stars to evolve, and the SFR reached values up to 10'' Mq yr"' . 

|Pipinoetti ( l201lh adopted models developed bv lCaluraet"al ( l2008h for elliptical galax- 
ies. The models comprise dust contribution from different stellar sources, a QSO wind 
and dust grain growth in the ISM. A model galaxy as massive as 10'^ Mq was applied to 
Jl 148+5251. Although the predicted SFR exceeds 3 x 10^ Mq yr"', the observed large dust 
masses could only be reproduced with a strong contribution from dust grain growth in the 
ISM in addition to dust produced by SNe and massive AGB stars. The results are strongly 
affected by a very high assumed dust destruction due to SN shocks in the ISM. Moreover, 
the models assume dust from SNe la, even though there is no clear evidence for significant 
dust production by these Types of SNe (see Sect. [533}. The dust contribution by the QSO 
wi nd of a few times 10^ Mp is in sufficient. 

iDwek and Cherchnefdfcoilh constructed scenarios comprising star-formation histories 
with a dominant dust production by either AGB stars or SNe. An average dust yield of about 
0.15 Mq (Cherchneff and Dwek, 2010) for all SNe, independent of their progenitor mass, 
was assumed. Only in cases of short-duration and intense bursts with SFRs in excess of 10** 
Mq yr"', SNe are found to be sufficient to produce the observed dust masses^ 

Pertaining to the discrepancies of the claims made in previous models, iGall et all ( 12011^ 
ascertained the impact of diverse astrophysical conditions governing the evolution of the to- 
tal dust mass in galaxies. The model takes into account AGB stars, different Types of core 
collapse SNe, and stellar yields from different groups. The dependence of stellar lifetimes 
on stellar masses and dust destruction due to SN shock interactions are considered. A simple 
treatment to estimate the impact of the formation of a supermassive black hole is introduced. 
It is shown that the amount of dust reached in galaxies and the significance of the contribu- 
tion by either AGB stars or SNe strongly depend on the assumed mass of the galaxy and is 
sensitive to the interplay between the IMF, the SFR, the dust production efficiency of SNe 
and the degree of dust destruction. Overall, larger dust masses are achieved with increasing 
mass of the galaxies or IMFs biased towards higher stellar masses. The calculations show 
that for increasing mass of the galaxies (and fixed SFR and IMF) either an increasing degree 
of dust destruction or a lower SN dust production efficiency can be accommodated. 

iGall et a il( l2011bl) identified plausible scenarios capable of reproducing the large ob- 
served dust masses for different QSOs at z > 6. They found that large quantities of dust can 
be generated in QSOs as early as 30-170 Myr after the onset of the starburst if the SFR of 
the starburst is > lO-' Mq yr"'. An initial gas mass of the galaxy of about 1-3 x 10" Mq was 
found to be sufficiently large. However, SNe are required to be very efficient while at these 
early epochs AGB stars contribute only marginally (see also this work. Sect. 17.21 FigjUl. 

It is worth stressing that the predictions of chemical evolution models partly reflect 
the initial assumptions made about physica l conditions and p rocesses. For example, the 
prediction of a major AGB star contribution jValiante et aii[2009() can be reproduced by the 
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models of lPwek and Cherchnefj bOllh and lOall et all ll2011allbl) when similar galaxy mass, 
SFR, IMF and dust production of stellar sources are assumed. 

Chemical dust evolution models face many uncertainties constituting various caveats 
in the aferomentioned approaches. Apart from the ambiguous dust production by stellar 
sources discussed in this review, the unknown formation and evoluti on of the QSOs and the 
assoc i ated star f o rmatio n history, pose a fundamental problem ("e.g.. lDwek and Cherchnefn , 
I2OIII : iGall et"5l2011ah . The composition and shape of dust grains produced by stellar 
sources or reprocessed in the ISM are relatively unknown. Commonly, carbon or silicate 
type dust are assumed, but other grain species might be considered as well (e.g., iDwekl 

Furthermore, destruction and growth processes of dust in the ISM, which decisively 
influence the hfetime of dus t grains, are poorly understood (e.g., lLiffman and Clavtonll 19891 : 
iMcKeel \l9S^ Ijonesl |2004 iDwek et all |2007[) . Estimates of grain lifeti mes resulting from 
calculations for the Milky Way range betwe en 100 and 1000 Myr (e.g.. Ijones et ai ll994h . 
although these predictions are uncertain (e.g.. |jones and Nuthll201lh . However, taking these 
lifetimes for granted, it has been s hown that AG B stars together with SNe can account for 
only ~ 10% of the interstellar dust llDraineLl2009h in the ISM of the Milky W ay and leads to a 
'miss ing dust source problem' in the LMC similar to that in high-z galaxies jMatsuura et all 
l2009h . 

Major outstanding qu estions for high-z QSOs thus include the effects of (i) dust grain 
_£rowth in the ISM (e.g., I PraineL l2009l : |Michalowski et al', '2010a'; 'Dwek and Cherchnef^, 
2011) . which is the preferred scenario in the Milky Way (e. g., Zhu k ovska et al, 2008; Draine, 
2009h. and (ii) d ust formation in QSO outflows (e.g.. lElvisetai I2OO2I ; IPipino et aj I20T1I : 



GaU etalLl201M . 



9 Summary and Conclusion 

This review has been devoted to dust formation by massive stars progenitors, including their 
role as dust producers in galaxies with emphasis on the early Universe. At very high redshift 
(z ^ 6) the minimum stellar mass of potential dust sources is ~3 Mo (see Fig.ffli. 

In Sect[2]we have discussed the many different channels in which stars with masses ~3 
Mo can evolve towards their dust producing end stages. Stars at the lower mass end (3-8 
Mo) evolve to AGB stars and release dust through intense mass-loss, most efficiently at the 
very end stages of evolution. Observationally, mass-loss rates can be obtained and linked to 
dust-mass-loss rates via gas-to-dust mass ratios, but it is difficult to determine the mass of 
the AGB stars. Thus, current information about the mass dependency of gas and dust yields 
from AGB stars relies on theoretical models (Sect|4ll. The current state of affairs indicates 
that 3^ Mo AGB stars are dominant dust producers among the AGB stars due to a higher 
dust production efficiency convolved with the IMF (see Fig. |7]l. However, the theoretical 
models are still rather simplistic, because e.g., the driving mechanism for mass-loss, the 
chemical composition of the atmospheres and hence the details of dust formation are poorly 
understood. 

Stars more massive than ~ 8 Mo explode as SNe but there is no clear one-to-one corre- 
spondence between the mass of the progenitor, the Type of the SN and the amount of dust 
produced. The most common Type of observed SNe with reported dust masses are Type II-P 
SNe. However, as discussed in Sec.|5]there are other promising channels (i.e., ECSN, Type 
Iln SNe, LBV stars) for producing significant amounts of dust. On the other hand, there is 
no strong indication or evidence for significant dust formation in other SN Types (i.e., la, lb. 
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Ic, and lib). Based on the currently available sample of observed SNe and SNRs (Tables[3] 
and|4ll there is a correlation between the observationally derived dust temperature and the 
amount of dust (see Fig.[2j3). The amount of inferred cold dust (< 100 K) is higher than that 
of hot dust (> 100 K). While there is no evidence for cold dust in SNe at early epochs, small 
amounts of hot dust are present in SNRs (see Fig.|2^). 

The results of theoretical models developed for dust formation in SNe are not in agree- 
ment with observations, predicting higher dust masses to be formed than observed. To ac- 
count for this discrepancy, models have been developed to investigate the effect of a SN 
reverse shock initiated by the collision of the SN forward shock with the ISM. However, the 
timescales for destruction (up to 10'* years) are too long to affect the dust masses in SNe, 
when observed at earlier epochs. 

In most models, the investigated SNe are considered to be Type II-P SN-like in the 
mass range 12^0 Mq, but the many complex physical and chemical processes involved 
are not yet well understood. Stars in the mass range 8-12 Mq have not been considered in 
the theoretical models for dust formation. The reasons are the lack of calculations of stellar 
yields due to the very complex evolution of these stars as discussed in Sect. l2.4l Substantial 
work has been devoted to model PISNe and zero-metallicity stars in the mass range 12^0 
Mq, even though so far there is little observational evidence for such stars. On the contrary , 
there is evidence for supersolar metallicity in very dusty high-z QSOs (e.g. jFan et ail2003l : 
iFreudlinget ai l2003l :l Juarez et al .20091 However, there are no models of dust production 
by SNe at supersolar metallicity. 

In Sect[6]we have parameterized the dust production efficiency of AGB stars and SNe 
using currently available observations and models. For operational purposes, we have ascer- 
tained a 'low', 'high' and 'maximum' case for the dust production efficiency of SNe. Using 
these efficiencies we have evaluated the total dust productivity of AGB stars and SNe in the 
mass interval 3^0 Mq for five different IMFs. We find that the dust production efficiency 
for AGB stars and SNe exhibit a decreasing tendency with increasing progenitor mass (see 
Fig. |4]i. The dust productivity of stars between 3^0 Mq is sensitive to the choice of IMF 
(Sect l7.2l i. This is more pronounced when stars between 8^0 Mq form dust with 'maxi- 
mum' SN efficiency. The contribution from AGB stars to the total dust productivity prevails 
when SNe are assumed to produce dust with a 'low' efficiency, but becomes insignificant 
for 'high' to 'maximum' SN dust production efficiency and for IMFs biased towards high 
stellar masses. The SN mass ranges of 8-12 Mq and 12-20 Mq are equally important and 
together dominate the dust production from all SNe between 8^0 Mq. Dust produced by 
stars between 20-40 Mq depends on the fractions of the various types of CCSNe. 

The present situation at high redshift is that large dust masses (> 10^ Mq) have been in- 
ferred from detection of thermal dust emission at sub-millimeter and millimeter wavelengths 
in the most distant QSOs. The issue about the origin of these high dust masses remains elu- 
sive and poses many unanswered questions. The most intriguing but also most uncertain 
ones concern the main dust sources, dust destruction and growth processes, the star for- 
mation history and the formation and evolution of the galaxies itself. On the other hand, 
deducing the dust masses from observations is challenging and afflicted with uncertainties 
(Sect.[8]l. Biases arising from assumptions about the IMF, the duration of the star formation, 
the absorption coefficient, the dust emissivity or the possible AGN heating could possibly 
lead to revisions of the dust masses currently quoted. Finally, the outcome of galactic chemi- 
cal evolution models also exhibit important differences, e.g., regarding the adequacy of AGB 
stars to produce sufficient amounts of dust at high redshift. The main reason for this is that 
the models greatly differ with respect to the assumptions made for the mass of the galaxy. 
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dust contribution from stellar sources as well as the treatment of the star formation history 
(Sect.IOl. 

Future progress on the theoretical side should involve more refined models for (i) dust 
formation in stellar outflows and (ii) dust evolution models for galaxies. For example, mod- 
els involving a more detailed account of the diverse possible end-stages of massive stars 
(e.g., ECSN, Iln, LBV) and the physical conditions leading to dust formation (SN ejecta, 
stellar winds in AGB stars, LBV outbursts) must be developed. Dust formation models 
should also allow for supersolar metallicity. Future galaxy models must be self-consistently 
connected to the mechanisms of the formation and evolution of the host galaxy impacting 
the evolution of the dust. Alternative dust sources complementing the stellar dust production 
such as dust grain growth in the ISM or dust production in QSO outflows might also be of 
relevance. 

From an observational point of view it is necessary to increase the currently small sample 
of observations of dust in SNe and extend it to different SN Types with indications of dust 
formation. Concerning the reliability of the observationally derived properties in the high-z 
QSOs, studies enlightening the possible influence of the uncertain parameters on e.g., the 
derived dust mass, SFR or stellar mass will be of relevance. 
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